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Two-dimensional  flow  of  an  incompressible,  inviscid  fluid  in  a  region 
with  a  horizontal  bottom  of  infinite  extent  and  a  free  upper  surface  is 
considered.  The  fluid  is  acted  on  by  gravity  and  has  a  non-diffusive, 
heterogeneous  density  which  may  be  discontinuous.  It  is  shown  that  the 
governing  equations  allow  both  periodic  and  single-crested  progressing  waves 
of  permanent  form,  the  analogues,  respectively,  of  the  classical  cnoidal  and 
solitary  waves.  These  waves  are  shown  to  be  critical  points  of  flow  related 
functionals  and  are  proved  to  exist  by  means  of  a  variational  principle. 
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SIGNIFICANCE  AND  EXPLANATION 


Th«  research  in  experimental  and  theoretical  hydrodynamics  in  the  last 
few  decades  has  indicated  that  solitary  waves  play  a  special  role  in  the 
evolution  of  general  disturbances  in  fluids.  Still,  the  investigation  of 
solitary  waves  and,  in  particular,  the  use  of  variational  principles 
associated  with  these  waves  is  far  from  complete.  While  variational 
principles  for  surface  waves  in  fluids  of  constant  density  have  been  discussed 
in  the  literature,  the  existence  proofs  given  here  appear  to  be  the  first 
rigorous  use  of  critical  point  theory  to  obtain  surface  waves.  Moreover,  we 
treat  a  class  of  density  profiles  not  heretofore  included  in  an  exact  theory. 

In  this  report  we  treat  a  two-dimensional  flow  of  an  incompressible , 
inviscid  fluid  in  a  region  with  a  horizontal  bottom  of  infinite  extent  and  a 
free  upper  surface.  The  fluid  is  acted  on  by  gravity  and  has  a  non-diffusive, 
variable  density  which  may  be  discontinuous.  It  is  shown  by  means  of  a 
variational  principle  that  the  governing  equations  allow  both  periodic  and 
single-crested  progressing  waves  of  permanent  form,  the  analogues, 
respectively,  of  the  classical  cnoidal  and  solitary  waves.  The  solitary  waves 
are  obtained  from  periodic  ones  as  the  periods  grow  unboundedly.  All  of  the 
waves  obtained  have  elevated  streamlines  and  have  speeds  greater  than  the 
critical  speed  associated  with  the  ambient  density.  Further,  the  amplitudes 
are  shown  to  be  exponentially  decreasing  away  from  the  crest. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
sumary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


A  variational  approach  to  surface  solitary  waves 


R.  B.  L.  Turner 

INTRODUCTION 

This  paper  is  concerned  with  two-dimensional  flow  of  an  incompressible,  inviscid  fluid 
in  a  region  with  a  horizontal  bottom  of  infinite  extent  and  a  free  upper  surface.  The 
fluid  is  acted  on  by  gravity  and  has  a  non-diffusive  variable  density  which  may  be 
discontinuous.  It  is  shown  that  the  governing  equations  allow  both  periodic  and  single- 
crested  progressing  waves  of  permanent  form,  the  analogues,  respectively,  of  the  classical 
cnoidal  and  solitary  waves.  Moreover,  solitary  waves  are  shown  to  arise  from  periodic  ones 
as  the  period  grows  unboundedly.  A  survey  of  earlier  work  on  steady  waves  in  stratified 
fluids  and  references  to  the  literature  are  given  in  [11  and  [2].  The  work  on  surface 
waves  in  fluids  of  constant  density  has  a  much  longer  history,  going  back  to  the  middle  of 
the  nineteenth  centuryi  see  [2],  [3],  [4]  for  references  and  accounts  of  the  development  of 
the  subject. 

The  problem  treated  here  is  close  to  that  examined  by  Ter-Krikorov  [5]  who  treats  a 
smoothly  varying  density,  decreasing  with  height,  and  allows  a  free  or  fixed  upper 
surface.  He  shows  that  from  each  vertical  mode  of  a  linearized  flow  problem  there  is 
bifurcation  to  a  wave  of  arbitrarily  prescribed  horizontal  period,  including  that  of 
"infinite  period",  i.e.,  a  solitary  wave.  The  methods  used  are  close  to  the  perturbation 
technique  of  Friedrichs  and  Hyers  [6]  who  gave  an  alternate  proof  to  that  of  Lavrentiev  [7] 
for  the  existence  of  small  amplitude  surface  solitary  waves.  The  techniques  used  here  are 
variational  in  nature.  They  are  an  outgrowth  of  the  work  of  Bona,  Bose,  and  Turner  [2]  on 
smoothly  stratified  flows  in  regions  with  fixed  upper  and  lower  boundaries  and  are 
particularly  close  to  the  methods  used  by  the  author  in  [1],  wherein  we  considered  two 
fixed  boundaries,  but  allowed  a  discontinuous  density.  The  present  paper  is  based  on  the 
observation  that  the  free  surface  can  be  treated  as  an  additional  discontinuity  at  which 
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the  density  drops  to  zero.  The  estimates  and  proofs  required  here  are  similar  to  those  in 
[1],  differing  mainly  in  the  behavior  near  the  free  surface.  To  shorten  the  length  of  the 
presentation  we  will  frequently  make  reference  to  proofs  in  [1],  pointing  cut  how  the  free 
surface  is  accomodated. 

Variational  principles  satisfied  by  free  surface  flows  of  constant  density  have  been 
described  by  a  number  of  authors  (cf.  e.g.  [8]-[13]  and  their  lists  of  references).  In 
[9],  (101  and  [13]  a  case  is  made  for  the  use  of  dynamically  invariant  quantities  in  a 
variational  characterisation  of  a  flow  as  a  step  toward  a  treatment  of  stability  in  the 
spirit  of  Liapunov.  Similar  ideas  have  been  carried  through  for  the  Korteveg  de  Vries  and 
other  model  equations  (of  [14],  [15],  [16]).  It  appears  difficult  to  base  a  rigorous  proof 
of  the  existence  of  steady  wave  solutions  of  the  Euler  equations  on  the  dynamic  principle 
given  in  [13].  Here  we  obtain  waves  using  a  differant  principle  and,  to  our  knowledge, 
ours  is  the  first  rigorous  use  of  a  variational  procedure  to  obtain  periodic  and  solitary 
surface  waves.  Garabedian  gives  a  critical  point  principle  for  periodic  surface  waves,  but 
it  appears  that  his  appeal  to  Morse  theory  needs  further  justification.  It  can  be  shown 
that  the  functional  he  uses  is  not  uniformly  positive  definite  at  the  origin,  as  claimed. 
The  principle  used  here  is  Lagrangian  in  character  and  reduces  to  a  constrained  variational 
problem.  The  separate  functionals  used  are  not  constants  of  time  dependent  motions  and  so 
the  method  does  not  immediately  suggest  a  means  for  establishing  stability.  However,  we 
feel  that  a  further  understanding  of  the  variational  structure  of  the  problem  will  be 
useful.  It  should  be  noted  that  in  a  related  problem,  that  of  vortex  flow,  principles 
allied  to  dynamics  have  been  successfully  used  to  establish  the  existence  of  steady  flows 
(cf.  [17]  and  references  given  there). 

Here  a  steady  wave  will  correspond  to  a  critical  point  of  a  "displacement"  functional 
on  a  manifold  of  prescribed  "kinetic  energy"  R.  For  the  classical  solitary  wave  what 
corresponds  here  to  R  almost  certainly  takes  valuss  in  a  finite  range  0  <  R  <  R  <  • 

(the  continue  found  in  [4]  have  this  property).  Thus  it  is  to  be  expected  that  with 
variable  density  the  range  of  R  will  be  finite.  This  limitation  is  reflected  in  a  lack 
of  coercivity  in  the  analytical  problem  derived  here.  Our  method  of  treatment  involves 


-2- 


Introducing  artificial  coarcivity  and  coabining  thia  with  estimates  on  tha  site  of  tha 
aolution  of  the  altered  problea.  The  outcome  for  the  original  problea  ia  a  reatriction 
0  <  R  <  R  on  allowable  energiea  with  the  aise  of  R  buried  aaong  elliptic  eatiaatea. 
While  it  ia  difficult  to  caopare  the  aize  of  tha  aolutiona  obtained  here  with  those  of  Ter- 
Krikorov,  they  auat  both  be  considered  to  have  small  solitude.  However,  the  variational 
approach,  with  iaproved  eatiaatea,  could  provide  finite  amplitude  waves.  Xpert  from 
treating  less  regular  densities  than  considered  in  [5]  we  can  show  that  in  the  presence  of 
a  free  surface  there  are  always  waves  of  elevation  with  aaplitudea  decreasing  away  from  the 
crest.  On  the  other  hand,  Ter-Krikorov'a  techniques  give  an  explicit  asyaptotic  fora  for 
the  wave  near  a  bifurcation  froa  a  parallel  flow,  though  not  uniformly  in  the  horizontal 
variable,  and  are  applicable  to  bifurcation  froa  higher  vertical  modes,  not  covered  by  the 
treatment  here.  It  should  be  noted  that  in  the  case  of  constant  density  it  has  been  shown 
that  faailiea  of  cnoidal  and  solitary  wavaa  exist  which  include  saall  amplitude  waves  and 
tha  Stokes  wave  with  a  sharp  crest  having  a  120*  opening  (cf.  [4],  [IB],  [19]. 

The  organization  of  the  present  paper  is  as  follows.  In  section  1  the  physical 
situation  is  described  and  the  relevant  mathematical  aquations  set  down.  The  analytical 
problems  are  posed  and  the  stain  results  are  described  in  theorem  1.2.  The  reaainder  of  the 
paper  is  devoted  to  establishing  these  results.  In  section  2  a  variational  formulation  is 
given  and  an  "extended*  problea  with  artificial  coarcivity  is  solved,  in  section  3 
estimates  for  the  solution  of  the  extended  problea  are  given.  Section  4  contains  estimates 
which  establish  that  solutiona  of  the  extended  problem,  when  restricted  to  have  small 
energy,  solve  the  original  flow  problea  in  the  periodic  case.  Section  5  deals  with 
exponential  decay  of  wave  amplitudes  away  froa  the  crest  and  tha  existence  of  solitary 
wavea  aa  the  limit  of  periodic  ones  when  the  periods  increase  indefinitely. 


1.  STEADY  PLOWS  WITH  A  PRES  SURFACE 


Here  we  briefly  describe  the  passage  fra*  a  physical  isodsl  of  wave  notion  to  a 
boundary  value  problem  for  a  partial  differential  equation.  For  a  more  coaplete  discussion 
we  refer  to  (1]  and  [2].  Consider  a  heterogeneous,  incompressible  fluid  acted  upon  by 
gravity  and  restrict  attention  to  flows  which  are  two-dimensional.  That  is,  assume  all 
quantities  depend  only  on  a  horisontal  coordinate,  a  vertical  coordinate,  in  the  direction 
of  gravity,  and  on  tine.  The  fluid  is  further  assuned  to  be  inviscid  and  non-diffusive, 
the  latter  property  to  be  elaborated  in  the  following  paragraphs.  While  our  interest  is  in 
wave  patterns  which  progress  horisontally  at  a  fixed  velocity  c,  we  can  remove  the  time 
dependence  by  considering  Cartesian  coordinates  referred  to  a  moving  crest  of  e  wave.  It 
is  then  possible  to  seek  a  steady  flow  in  a  region  which  is  independent  of  time.  The 
region  will  have  to  be  determined  as  part  of  the  solution  of  the  problem.  However,  we  do 
assume  the  flow  is  over  an  infinite  horisontal  bottom.  Let  x  be  a  coordinate  in  the 
horisontal  boundary  and  y,  a  vertical  coordinate  chosen  so  that  the  bottom  boundary  is 
at  y  »  -1  and  so  that  the  acceleration  due  to  gravity  is  represented  by  (o,-g)  with 
g  >  0.  The  fluid  ia  assumed  to  have  a  free  surface 

y  -  h(x)  >  -1  (1.1) 

which  is  to  be  determined.  To  begin  with  we  require  h  to  be  continuous  and  satisfy 

lim  h(x)  -  0  (1.2) 

The  fluid  is  then  assumed  to  occupy  the  region 

Fft  -  {(x,y)|-1  <  y  <  h(x),-w  <  x  <  »}  (1.3) 

For  (x,y)  «  F.  let 
n 

a  -  (o.v)  (i.4) 

where  U(x,y)  and  V(x,y),  respectively,  are  the  horisontal  and  vertical  components  of 
fluid  velocity  in  a  moving  frame.  A  flow  is  sought  which  is  steady  in  the  moving  frame  so 
the  Euler  equations  take  the  fora 

P<S  •  v>a  -  -Vp  -  gpgj  in  Fh  (1.5) 

where  p  “  p(x,y)  is  the  density,  p  is  the  pressure,  and  e_  -  (0,1).  The  condition  of 


incompressibility  becomes 


div  g  ■  0  (1.6) 

and  the  condition  of  nondif fuaivity  antaila 

g  •  Vp  -  0  (1.7) 

throughout  Fh.  Supplementing  aquations  (1.5),  (1.6)t  and  (1.7)  we  have  the  following 
boundary  conditions.  Tha  fluid  should  not  panatrata  tha  horizontal  bottom;  i.a., 

V(x,-1)  -  0,  —<*<«.  (1.8) 

On  tha  fraa  surfaca  S  «  {(x,y)jy  ■  h(x) <  x  <  "}  there  ara  two  conditions!  a 
kinaauitic  condition  requiring  the  velocity  to  be  tangential  and  a  second  requiring  tha 
pressure  to  be  zero;  i.a. 


-h  V  +  V  *  0  on  S 
x 


(1.9) 


and 


p  ~  0  on  S  .  (1.10) 

Hera  hx  denotes  a  derivative.  Finally,  conditions  oust  be  specified  at  x  -  ±«.  We  can 
specify  a  velocity  distribution  at  infinity  and  here  take  the  simplest  case  of  a  wave 
propagating  in  fluid  which,  in  “laboratory"  coordinates,  is  at  rest  at  infinity.  For 
coordinates  based  in  a  wave  moving  to  the  left  with  velocity  c  the  condition  becomes 


lim  (0,V)  -  (c,0) 
x*t« 


(1.11) 


Further,  the  density  in  the  “undisturbed  region*  is  specified  and  we  take  it  to  depend  only 
on  y.  Thus 


lim  p(x,y)  •  p(y) 
x*±« 


(1.12) 


where  p(y)  is  a  decreasing  function  of  y  for  -1  <  y  <  0,  normalized  to  satisfy 
P(0)  -  1. 

The  problem  that  will  occupy  us  is  to  find  c,  h,  g,  p  and  p  satisfying  (1.1), 
(1.2),  (1.5)— (1.12).  The  regularity  required  of  the  functions  will  be  made  more  precise  in 
what  follows.  For  now  we  continue  the  discussion  at  a  formal  level.  In  particular,  we 


••sum  to  begin  with  that  p  ia  continuously  differentiable  although  we  ultimately  want  to 
allow  discontinuous  densities. 

The  conditions  (1.6)  and  (1.7)  imply  that  there  is  a  pseudo- stream- function 
♦  -  Vtx.y)  such  that 


14 

3y 


p1/2u, 


-  ¥  ■  PV2* 

9x 


(1.13) 


and  that  p  is  a  function  of  V,  p(f).  From  Bernoulli's  theorem  the  total  head  is 
constant  along  streamlines,  hence  on  level  sets  of  f.  Thus 

H  -  p  +  i  P(D2  ♦  v2)  +  pgy  -  H(  *)  (1.14) 

Eliminating  p,  using  the  two  components  of  the  vector  equations  (1.5).  one  finds  that  4> 
satisfies  the  equation 

A*(x,y)  +  gy  jj  -  fj  (1.15) 

(cf.  Dubreil- Jacotin  (20],  Long  (21]  and  Yih  [22]).  If  the  density  p  is  specified  as  a 
function  of  and  if  the  dynamics  are  specified  by  giving  H(4>),  then  (1.15)  is  a 
semi linear  elliptic  equation  for  f(x,y).  Any  solution  of  (1.15)  gives  rise  to  a  solution 
of  (1.5)  with  U,  V  and  p  obtained  froa  (1.13)  and  (1.14).  For  now  we  leave  aside  the 
question  of  boundary  conditions  on  the  top  and  bottoa  of  Fj,  and  examine  the  implications 
of  (1.11)  and  (1.12). 


In  the  search  for  a  solitary  wave,  a  disturbance  which  should  be  of  essentially  finite 
extent,  it  is  natural  to  ask  that  for  large  x,  t(x,y)  should  approach  a  pseudo-stream- 
function  corresponding  to  a  flow  with  velocity  (c,0)  in  a  stream  of  density  p(y).  Thus, 
letting 


we  require 


♦<y> 


c  f  P1/2(.)ds 


(1.16) 


♦(x.y)  -  Wy>  ♦  0  as  |x|  ♦  »  ,  (1.17) 

The  density  p  (we  use  the  same  symbol)  associated  with  the  stream  coordinate  ♦, 
coapatible  with  that  already  specified  at  x  *  ±mt  is 

P(»)  -  P(y(*)>  (1.18) 

where  y ( ♦)  is  the  function  Inverse  to  {(y).  The  expression  (1.18)  is  taken  as  the 
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definition  of  p(f).  Similarly,  the  total  head  associated  with  the  stream  coordinate  t 
is  taken  to  be 


H<*)  -  pty< i  p(*)c2  ♦  p(*)gj(*) 


where  p  is  hydrostatic t 


P 


p ( s ) da  . 


(1.19) 


(1.20 


With  p(1>)  and  H(f)  given  we  can  now  seek  a  solution  t(x,y)  of  (1.15).  Mote  that  for 

O 

any  constant  c,  $(y)  (cf.  1.16)  is  a  soltuion  of  (14).  We  call  it  a  trivial  solution. 

The  flows  examined  here  will  have  free  surfaces  and  discontinuities  in  velocity 

occurring  along  certain  streamlines.  in  order  to  deal  with  a  problem  in  a  fixed  domain  and 

to  confine  irregular  behavior  to  coordinate  lines  we  replace  (1.15)  by  an  equation  for  y 

in  semi-Lagrangian  independent  coordinates  x  and  f.  The  interest  here  is  in  flows 

without  reversal*  i.e.  with  0  >  0  or  >  0  and  so  an  inversion  yielding  y(  x,  Ji)  is 

oy 

reasonable.  Corresponding  to  (1.15)  is  the  equation 


(1.21) 


(cf.  [1],  [5]),  obtainable  using  the  relations 

*x  +  Vx  ’  °'  V*  “  1  <1,22> 

O  O 

Naturally,  y  -  y(*),  the  inverse  of  t(y),  is  a  solution  of  (1.21)  and  we  refer  to  it  as 
trivial.  Corresponding  to  (1.17)  is  a  condition 

y(x,f)  -  y(f)  ♦  0  as  |x|  ♦  -  ,  (1.23) 

which  is  imposed  uniformly  in  ♦  .  <  ♦  <  0  where 


♦(-1) 


/  pV2(s)ds 
0 


(1.24) 


The  condition 


y(x,<i_1)  -  -1  (1.25) 

replaces  (1.9),  and  from  (1.22)  it  is  clear  that  (1.9)  will  be  satisfied  merely  by  choosing 
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I 


h(x)  -  y(x,0)  . 


2  2  2  -2 

Since  p(U  +  V  )  -  (1  +  y‘ )y^  ,  the  pressure  condition  at  the  fluid  surface  become 


p(x,0 )  «  constant 


(1.27) 


1  '  r 

p(x,*)  -  H(*)  -  r:  ( - r-5)  -  p(*)gy(x,*)  , 


using  (1.14).  Suppose  that  a  continuously  differentiable  p(y)  is  given  and  that  for  any 


fixed  speed  c,  pit)  and  H(40  are  defined  by  (1.18)  and  (1.19).  It  is  then  possible  to 


interpret  the  foregoing  equations  and  conditions  in  a  classical  sense.  For  a  density  p 


with  possible  discontinuities  a  meaning  must  be  given  to  equation  (1.21)  and  we  do  that 


Let  c  '  on  a  domain  consist  of  functins  with  continuous  derivatives  through  order 


k,  each  satisfying  a  Httlder  condition  with  exponent  8.  For  a  more  complete  description 


of  spaces  see  Ml,  section  1.  Suppose  p(y),  given  initially,  satisfies 


i)  p  is  nonincreasing  on  1-1,0] 


ii)  pec*  for  some  8  >  0  except  for  jumps 


(1.29) 


at  points  n1  i  -1  <  n1  <  n2  <  ...  <  nH  <  0 


where  p  is  continuous  from  the  right. 


The  corresponding  function  p(l>)  from  (1.18)  will  have  discontinuities  at  points 


1  <■  1  10  N+1 

if  “  ♦(!)  ),  j  -  1,2, ...,N.  Extending  the  domain  of  f  to  n  "  -1  and  n  “0  (so 


that  »  V(n°)  -  )  we  let 


O’1  -  R  *  (*i,Vi+1)»  j  -  0,1,.  ...N 


and  D  »  *  x  for  0  <  a  <  1  let 


Ha  -  {y|y  e  C0' 1  (D)  ,y  e  C°(D),y  €  C1,a(Dj), 


and  y^  >  0  in  D  ,  j  ■  1,2,...,N> 


Since  (1.21)  has  a  divergence  form,  the  notion  of  a  weak  solution  where  p  and  H  are 


smooth  is  described  in  a  standard  way  (cf.  [1]  (1.13)).  A  weak  solution  can  also  be 


defined  with  the  presence  of  discontinuities  in  p  and  H  and  is  equivalent  to  the 
following  definition  which  stresses  a  physical  aspect  of  the  flow. 

Definition  1.1.  He  call  y  e  M*  a  solution  of  (1.21)  if  and  only  if 
i)  y  is  a  weak  solution  of  (1.21)  in  Dj  for  j  -  1,2,...,N. 

(ii)  The  pressure  p(x,4>)  computed  from  (1.2S)  in  each  D^,  j  “  0,1, 2,, ..,N  is  the 
restriction  to  D j  of  a  continuous  function  p  defined  on  D. 

The  solitary  wave  problem  can  be  posed  as  follows. 

Problem  P^.  Find  a  c  >  0  and,  for  some  a  >  0,  a  nontrivial  (1.32) 

function  y  e  m“,  satisfying  (1.21),  (1.23),  (1.25)  and  (1.27). 

In  the  course  of  solving  the  solitary  wave  problem  we  will  solve  the  corresponding 

problem  for  periodic  waves. 

Problem  Pfc.  Find  a  c  >  0  and  a  nontrivial  y  e  M°,  y  a  (1.33) 

2k  periodic  function  of  x,  satisfying  (1.21),  (1.25)  and  (1.27). 

The  interpretation  of  a  solution  of  problem  P^  in  terms  of  a  solitary  wave  hinges  on 
the  condition  (1.23),  requiring  the  flow  to  be  "trivial"  at  x  -  *».  while  the  physical 
interpretation  of  a  solution  of  problem  Pk,  for  finite  k,  is  tenuous,  we  will 
nevertheless  find  solutions  which  are  2k  periodic  in  x  and  show  that  they  are 
exponentially  decaying  from  crest  to  trough.  Thus  they  might  be  viewed  as  a  train  of  waves 
with  quiescent  zones  of  approximately  trivial  flow  between  the  crests.  In  any  case,  the 
periodic  formulation  is  a  convenient  analytical  device. 

To  analyze  the  problems  P^  for  k  <  »  and  to  expand  the  class  of  density  profiles, 
we  associate  with  any  density  function  p  described  by  (1.29)  a  family  of  smooth  densities 
as  follows.  For  a  4g  which  is  positive  but  smaller  than  -  n^”1  for  j  ■  1,2,.,.,n 


and  for  each  6,  0  <  6  <  6  ,  1st  p.(y)  b«  a  nonincreasing  function  of  cists  C 


"  -t 

P»(y)  ■  Ply)  for  y  4  u  (hJ  -  4,nJ)  (1.34) 

J-1 

snd  sat  pQ (y )  -  p(y),  Tha  mannar  In  which  tha  jump  at  is  approximated  by  a  smooth 

transition  is  immaterial  to  our  estimates  and  so  we  need  not  specify  p ^  on  the  intervals 
( -  5,r^),  j  »  1,2, ...,N.  The  estimates  in  the  subsequent  sections  give  a  uniform 
picture  of  the  transition  from  smooth  to  discontinuous  densities. 

In  order  to  summarize  the  Min  results  of  the  paper  in  a  convenient  form  it  is 
necessary  to  anticipate  a  transformation  to  be  introduced  in  section  2  and  give  some 
additional  notation.  A  rescaled  stream  variable  is  given  by 

•»«♦»  .  j  —-SS- -  (1.35) 

0  cp  '■‘(y(s)) 

(cf.  2.14)  so  that  at  X  •  i“  the  streamline  with  label  n  has  height  n.  The  expression 

w(x,n)  -  y( x, ♦( n) )  -  n 

(cf.  2.2)  represents,  for  each  x,  the  displacement  from  its  position  In  a  trivial  flow, 
of  the  streamline  which  has  height  n  at  x  »  i*  (under  the  condition  1.23).  The 


equation  satisfied  by  w,  with  x  and  n  as  independent  coordinates,  is  given  by  (2.12) 
(where  x  •  x1  and  n  «  Xj).  Its  formal  linearisation  about  w  »  0  is 

f-  [|r  P(b)  ~  —  P(n)  ~j  -  g-^  s  in  **<-i,o)] 


z  »  0  at  n  ■  -1  i 


2  3z 

c  •z-  "  gs  at  ti  •  0 
on 


(cf.  4.22),  the  last  condition  arising  (as  we  shall  see  in  proposition  2.1)  from  the 

Bernoulli  condition  (1.27).  The  Sturm-Liouville  problem  obtained  from  (1.36)  by  letting 

2 

z  »  z(n)  has  a  least  eigenvalue  1  “  g/c  .  The  corresponding  value  of  c,  denoted  cs 
to  indicate  its  dependence  on  S  through  p  «  p^,  is  the  largest  speed  (correspondingly, 
the  lowest  spectral  point)  associated  with  (1.36).  This  speed  is  referred  to  in  the 
hydrodynamics  literature  as  the  "critical*  speed  or  the  speed  with  which  infinitesimal  long 
waves  travel.  Its  relevance  here  is  that  the  waves  we  obtain  have  propagation  speeds  which 


w-  ■••V 


art  "supercritical",  i.e.,  they  are  larger  than  c A  full  discussion  of  these  ideas  has 
been  given  by  Benjamin  [23]. 

The  main  results  of  the  paper  are  contained  in  the  following  theorem.  Its  proof,  and 
explicit  estimates  related  to  the  waves  obtained,  are  spread  over  the  ensuing  sections.  In 
particular  we  refer  to  propositions  2.1,  2.2,  4.4,  4.6,  4.8  and  S.2  together  with  remarks 
2.3,  4.9  and  5.3.  Note  that  a  discontinuous  density  is  allowed  when  the  parameter  S  “  0, 
while  a  solitary  wave  corresponds  to  the  case  k  «  ",  by  our  convention. 


Theorem  1.2.  There  are  positive  numbers  R,  x(R),  and  a  such  that  for  0  <  S  <  6Q, 

0  <  R  <  R,  and  k(R)  <  k  <  *  the  problem  for  p^  has  a  nontrivial  solution  y  in 

a 

M  corresponding  to  a  speed  c  with  the  following  properties 

1)  c  >  c4(1  -  C.jR4'^3)  1/>2  i.e.  c  is  "supercritical" 

2)  y  has  period  2k  in  x  (for  k  <  «•). 

The  streamline  displacement  w  satisfies 


3) 

0  k 
1  / 

p(ni  4^ —  dxdn  * 

R2 

-1  -k 

n 

4) 

w(x, n) 

>0  for  -1  <  n  < 

0. 

5) 

v(x, n> 

-  w(-x,n)  and  for 

0 

w(x,n)  >  w(x' ,n). 

6)  |w|  <  C^expl-Bx]  and  |Vw|  «  C^expl-Bx]  on  0  <  x  <  k  for  a  B  >  0. 

The  quantities  R,  k,  a,  Cy  C2  and  Cj  depend  on  p  in  (1.29);  C2  and  C3  also 
depend  on  R  and  B. 
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i 


2.  A  VARIATION  FORMULATION 


In  this  ssction  it  is  assumed  that  tha  density  is  smooth  and  given  by  (1.34)  for 
6  >  0,  but  the  subscript  is  often  omitted.  The  equation  (1.21)  is  formally  the  Euler 
equation  for  the  functional 


4<y) 


1  ♦  y  2 

//  k  (-7-*)  ♦  9P'(*)  ~  H*(*)y]dxd*  , 

d  Yii 


(2.1) 


where  primes  denote  derivates.  That  is,  the  condition  that  4  has  a  critical  point  (aero 
derivative)  at  the  function  y  is  expressed  by(1.21).  This  will  be  made  precise  at  a 
later  stage.  Note  that  the  first  term  in  the  integrand  in  (2.1)  is  merely  the  Dlrichlet 

O 

integral  in  the  new  coordinates.  As  noted  in  section  1,  y  “  y(  4)  is  a  solution  of 

(1.21) i  thus  it  is  formally  a  critical  point  of  4.  if  y(x,4)  is  another  solution  of 

(1.21)  then 

w(x,n)  *  y(x,J(n))  -  n  (2.2) 

is  formally  a  critical  point  of 


G(w) 


//  [5  c2p(n> 


■LM2-  » 

1  +wn 


gp'(n)  ~]dxdn 


(2.3) 


where 

a  »  {(x,n)l-«  <  x  <  •,  -1  <  n  <  0}  (2.4) 

(cf.  [1],  $2).  Here  n  is  the  stream  coordinate  introduced  in  (1.35)  and  w,  the 
vertical  itreamline  displacement,  (re  emphasize  that  for  a  nontrivial  flow  w  *  0,  in 
general*  i.e.,  only  at  x  -  ±-  does  one  require  that  the  streamline  with  label  n  have 
height  n.  Note  that  with  the  new  scaling  p{ n)  is  the  same  function  introduced  at  the 
outset,  describing  the  density  as  a  function  of  height. 

JUst  below  the  free  upper  surface  of  the  flow  the  density  has  the  positive  value  1. 
while  the  usual  model  for  surface  waves  implicitly  takes  p  -  0  in  the  atmosphere, 
retaining  only  the  Bernoulli  pressure  condition,  it  appears  necessary  for  a  workable 
variational  principle  in  the  present  context  to  explicitly  incorporate  a  drop  in  density  at 
tha  upper  surface.  Starting  with  p  (or  pg)  define 


„  p(n)»  -i  <  n  < 

p(n)  -  ■ 

o  »  n  «  o 


The  change  in  p  at  the  one  point  n  “  0  will  not  alter  the  firat  term  in  the  integrand 
in  (2.3)/  but  will  alter  the  aecond  tana.  Define 


Equivalently . 


2  2 

g(w)  -  //  cp(n)  j^‘-  ♦  gp'(n)  j-Jdxdn  . 


m  2 

G(w)  -  G(w)  -  /  g  — dx 


We  turn  now  to  the  periodic  version  of  the  wave  problee.  Let 

\  “  {(x,n)  I  |x|  <  k,  -1  <  n  <  0) 


G(w)  -  /  c2p(fl)  +  gp'(n) 


3 f 

It  will  often  be  convenient  to  use  the  notation  x,  “  x,  x  ■  n»  P,  -  —  and  f,  *  -r— 

*  2  1  1  ojp^ 

for  1-1/2  where 


pZ  ♦  p2 

..  .  1  11 _ 12 

f<P1,P2)  2  1  +  p„ 


(2.10) 


Still  proceeding  formally,  one  verifies  that  if  w  (assisted  to  vanish  when  x2  »  -1 )  is  a 

2 

critical  point  of  G,  then  with  X  m  g/c  , 

/  p(x  )f  (Vw)*s  »  X  /  p'(x  )w*  (2.11) 

ft  ^  1  Xi  A  * 


for  ellowable  "variations”  s.  Mere  and  in  the  sequel  a  repeated  index  is  understood  to  be 
ausaied  over  {1,2}.  The  equation  (2.11)  is  merely  a  weak  form  of 


^  ptx^f^Tw)  ■  Xp*(x2)w  in 


w  ■  0 


p(x2>f2<Vw)  ■  Xw 


•t  x2  -  -1 


at  x2  *  0  • 


(2.12) 


Ma  interrupt  tha  reformulations  of  the  problem  at  thia  point  to  show  that  a  nonsero  smooth 
solution  of  (2.12)  gives  rise  to  a  nontrivial  solution  of  problem  P^. 


2  i  - 

Proposition  2,1.  suppose  w  e  C  (0)  n  c  (Q).  If  w,  with  w^  >  -1,  satisfies  (2.12) 
for  some  X  >  0,  then  (cf.  2.2)  y  defined  by 

y(x,i>)  -  n(*)  ♦  w(x,n(*))  (2.13) 


where 


and 


n<*> 


* 

o 


da 

cp^tyls)) 


(2.14) 


c  -  g/A  (2.15) 

satisfies  equation  (1.21)  and  the  conditions  (1.25)  and  (1.27). 

Proof.  The  correspondence  between  the  elliptic  equation  (2.12)  for  w  in  Q  and  the 
equation  (1.21)  for  y  in  D  is  shown  exactly  as  In  section  2  of  (1).  also  the  assertion 
that  y  takes  the  value  -1  when  ♦  ■  ♦_1  follows  trivially.  To  see  that  the  pressure 
condition  (1.27)  is  satisfied  start  with 


f2(Pw)  -  Xw 

at  Xj  “  0  (note  p(0)  -  1)  and  express  the  relation  in  the  coordinates  x,  n> 


_  .  1  2  12 

(1  ♦  w_)2  c2 


(2.16) 


for  n  •  0.  To  compute  tha  pressure  from  (1.26)  the  value  of  H  for  n  ■  0  is  needed. 
Prom  (1.1*) 
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H|^.0  "  p(0)  +|c2  ♦  gy(0) 


Fro*  the  relations  (1.16)  and  (2.2),  evaluated  where  n  «  0  (or  <>  •  0)  and  p(0)  -  1, 

wx  "  *x 

and 


wn  -  V*n  -  1 


y^c  -  1 


Hence  the  express ion  (1.28)  for  tha  praaaure  becostss 


p(x,0)  ■  j  cJ  •  j  - -  c2  -  gw 

(1  +  V 

-  0 

whan  (2.16)  la  uaad,  completing  tha  proof  of  tha  proposition. 

To  obtain  nontrivial  solutions  of  Problsa  Pfc(k  <  »)  Tor  a  density  which  is  smooth  on 
“1  <  n  <  0  it  will  suffica  to  obtain  nonsaro  periodic  solutions  of  (2.12).  The  case  of 
discontinuities  in  density  on  -1  <  T>  <  0  and  the  resolution  of  problasi  for  solitary 

waves  will  be  handled  through  Halting  procedure*. 

The  equation  (2.11)  can  be  expressed  as 

F'(w)  -  XB'(w)  (2.17) 

where 


F(w)  •  /  p(x,)f (Vw) 


B(w) 


/ 

0, 


P’(x2) 


2 

w 

2 


(2.18) 


and  the  prises  in  (2.17)  denote  derivatives  (still  to  be  defined  in  a  suitable  space).  A 
tempting  approach  to  solving  aquation  (2.17)  is  to  consider 
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•up  B(w) 
F(w)-const. 


and  obtain  X  as  a  Lagrange  multiplier  associated  with  an  axtreaal  w.  However,  B  is 
unbounded  on  level  sets  of  F.  To  see  this  consider  the  case  that  p  S  1  and  let 

_  i  +  n 
/|*l  ♦  lnl1/2 

on  Q^.  Clearly  w(x,0)  f  L2i  however,  a  simple  estimate  shows  F(w)  <  ».  Nevertheless, 
as  shown  in  [1]  for  a  similar  problem  without  a  free  surface,  the  variational  approach  can 
be  salvaged  by  altering  f  in  the  integrand  of  F  where  1 9w|  >  r,  for  some  positive 
r,  and  then  showing  that  a  solution  obtained  by  a  variational  procedure  and  having  a 
suitably  restricted  "energy"  satisfies  |7w|  <  r.  The  approach  here  is  similar,  but  we 
also  taka  account  of  the  free  upper  surface  by  an  altered  functional. 

He  now  proceed  to  define  a  substitute  for  equation  (2.1?)  as  a  step  to  obtaining  the 
results  in  theorem  1.2.  He  refer  the  reader  to  section  2  of  [1]  for  proofs  of  some  of  the 

assertions  to  be  made.  He  are  still  considering  a  family  of  densities  p^,  0  <  4  <  SQ, 

(cf.  (1.34))  which  are  smooth  on  -1  <  n  <  0  and  which  could  reflect  a  rapid  change  in 
density  at  certain  levels  in  a  fluid.  A  similar  "smoothing"  for  p  near  n  »  0  will  be 

useful,  purely  as  an  analytical  device.  Let  t  -  x ( n )  be  a  nonincreasing  C2  function  on 

(-<•,0)  which  is  sero  for  x2  <  -1  and  which  satisfies  T(0)  ■  -1.  For  each 
e,  0  <  e  <  60  let  te(n)  “  t(n/c)  and  define 

Pj  e<n)  -  Pjtn)  +  Te(n)  (2.19) 

o  o  O 

so  that  p.(n)  -  lim  p,  .  As  earlier,  the  manner  in  which  p.  decreases  from  p.(-e) 

6  e>Q  o,e  S, t  6 

O 

to  0  on  [— c,0]  is  not  important,  though  we  do  require  that  pj  £  <  pj.  Now  extend  p^ 

to  0  <  n  <  1  as  an  even  function  and  p,  to  0  <  n  <  1  as  an  odd  function,  retaining 

"  o,c 

the  same  symbols  for  the  extended  functions.  The  extended  functions  will  be  in  class 
IB 

C  if  pj(0)  ■  Oi  otherwise  pj  is  merely  lApschits  continuous  at  0.  He  also 
require  the  function  s  defined  by 

«<n)  -  -sgn(n)  (2.20) 
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where  tha  signum  of  0  la  0.  Now,  suppose  {  -  C(t)  la  a  smooth  decreasing  cutoff 
function,  aqual  to  1  for  0  <  t  <  1  and  equal  to  0  for  t  >  2.  Raplacing 
in  (2.18)  will  ba 


»(x2,p1,P2) 


2  2 
.Pi  *  P, 


2  2 
Pi  ♦  P 


2  2 
Pi  ♦  P5 


>2(1 


*1  /r1  K2, 

.  *  "  -  v(— s— ) 


(2.21) 


whara  5  <p1#P2)  -  — -)  ao  that  a  «  (p2  +  p2)/2  whan  p2  ♦  p2  >  2r2.  Let 

^  2 

ai  -  3a/3pi  and  a^  -  3  a/Sp^p^.  Tha  function  a  ia  globally  convax  in  (p^.pj). 
uniformly  in  Xj,  for  r  aufficiantly  aaall,  and  aatiafiaa  tha  following  inaqualitiaa 
(provad  exactly  aa  in  [1],  leona  2.1). 


^(pj  p2)  <  a(p1,p2)  <  i  «2(p2  +  p2) 

°3<P*  ♦  <  •,?,  +  «2P2  *•  °y(P?  +  Pj> 

*1  +  *2  *  05(“lP1  *  W 

v(C?  ♦  C*>  < 


(2.22) 


Hara  oi  for  1  <  i  <  5  and  v  ara  poaitiva  conatanta  lndapandent  of  x2>  e,  pf  and 
P2- 

Lat 

&k  -  {(x^Xj)  |  |x1)  <  k,  |x2l  <  1)  (2.23) 

and  aat  Qm  ”  Q.  Lat  denote  tha  c"  functions  on  U  which  are  2k  periodic  in 

x,  and  have  support  not  containing  pointa  where  x2  -  *1.  Tha  symbol  denotea  tha 
elements  of  which  ara  even  with  reapect  to  x^  and  Xj.  Since  the  functiona  in  Ck 

and  c£  vanish  whan  x2  «  -1 ,  tha  Poincart  inequality  guarantees  that  the  expression 

lwlk  -  (/  |7w|V/2  (2.24) 

”k 

provides  a  norm,  and  tha  completions  of  tha  respective  spaces  in  the  norm  ara  denoted  by 
■  H^IQ)  and  H*  ”  H*(Q).  The  symbols  and  hJJ  will  also  ba  used  in  later 


•act Iona  to  danota  tha  raatrictlona  to  0  of  functions  in  thoaa  apacaa.  Thus  a  function 


in  H®(Q)  is  avan  in  x,.  For  w  6  H®  (Q)  dafina  tha  functionals 

A(w)  -  /  p^(x2>a(x2,7«)  (2.25) 

and 

B(w)  -  -/  p$  (Xj)  (2.26) 

where  integration  ia  over  f)^  hare  and  in  tha  raaaindar  of  this  section.  Tha  notation 
A* ,  B*  will  be  used  for  the  Frechet  derivatives  of  the  functionals  in  H®(  Q) . 

Our  subsequent  program  is  briefly  described  as  follows.  Me  show  that  for  each  R  >  0 

tha  aquation  (2.27)  below  has  a  solution  ( X,w)  with  w  »  w^  ^  £  positive  in  0^  and 
*  2 

normal lead  by  A(w)  •  2R  .  Restricting  attention  to  the  lower  region  B^  we  show  that  for 
e  converging  to  zero  through  a  suitable  subsequance  we  obtain  a  function  w  which  is 
smooth  on  and  satisfies  (3.23) ,  essentially  the  restriction  of  (2.27)  to  for 

e  «  0.  For  R  suitably  restricted,  |Vw|  <  r  and  we  can  show  that  ( A,w)  satisfies  the 
original  equation  (2.12).  All  estimates  obtained  are  uniform  in  6  >  0  and  k  >  0. 

Taking  limits  of  solutions  as  these  parameters  vary  we  obtain  the  desired  wave  forms. 

Until  we  consider  limits  involving  the  parameters  k,  6,  and  e,  we  selectively  suppress 
them.  Xt  is  assumed  that  k  <  «  until  Remark  5.3.  Through  the  penultimate  paragraph  of 
section  4  it  is  assumed  that  4  >  0,  while  for  the  remainder  of  this  section,  it  in  also 
assumed  that  e  >  0. 


Proposition  2.2.  For  each  R  >  0  the  problem 

A'(w)  -  \B*(w)  (2.27) 

0  A  2  <* 

has  a  solution  (X,w)  satisfying  1  >  0,  we  H^tB),  A(w)  -  2R  ,  and  w  >  0  in  B.  The 
function  w  is  characterized  by 


B(w) 


sup#  B( v+ ) 


vCTk 

A(v)-2R2 


where  v+  »  max(0,v). 

Proof.  The  proof  proceeds  as  that  of  Theorem  2.1  of  (11  with 
variational  procedure  leads  to  M($)  "0  for  all  +  e  H*(fl^) 


(2.28) 


H®  replacing  K^.  The 
and  for  a  suitable  X  >  0, 
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>:,•*** 


where  N  la  defined  fey 


M(«) 


/  tP<*2>»i(x2.VW)  {*- 


(2.29) 


The  derivatives  of  a  are  easily  seen  to  satisfy  (x2'”pi »P2*  "  ”*1 **2'p1 tP2*  *nd 
a1<"x2'pl'”p2)  ”  *1(x2'p1  'p2)'  whilB  aj  hM  th*  opposite  parity*  i.e.,  even  in  p, 
and  odd  in  <x2,p2).  further  p  and  S’  are  even  in  x2.  As  a  consequence  one  finds 
that  M  annihilates  all  functions  4  which  are  odd  in  or  in  x2  or  in  both.  Then 

M  is  zero  on  all  test  functions  and  the  remainder  of  the  proof  follows  as  before. 


Remark  2.3.  According  to  (24],  Theorem  6.3,  the  smoothness  of  w,  restricted  to  the 

1  B 

original  region  8,  is  limited  only  by  the  sstoothness  of  p.  Thus  with  p.  e  c  ,  w  is 

0 

2,6 

of  class  C  In  sach  subdomain  of  Q. 


Lsnwa  2.4.  The  multiplier  X  occurring  in  proposition  2.2  satisfies 

<  X  <  C2  (2.30) 

where  C, ,  C2  depend  on  the  total  variation  of  p( n) . 

Proof.  Since  p(0)  »  1,  it  follows  from  (2.27)  and  (2.22)  that 


i  -  4Llw)<J?  >  c  -Uadi  (2.31) 

<B’(w),w>  -/  p’w2 

Since  -/  p*w2  ■  2  /  pw  <  p(-1)  /  w2  (y^)2  <  CP(-1)  /  |7w|2  a  lower  bound  depending 
on  p(-1)  results. 

In  a  similar  fashion 


X 


<A* (w) ,w> 
<B’(w),W> 


20^0^  ^  A(w) 

bTw) 


“Si'*2 

5Twj 


(2.32) 


follows  from  (2.22),  (2.27),  and  the  charactarication  of  w.  The  last  quotient  in  (2.32) 

e  a  o 

can  only  become  larger  if  w  is  replaced  by  any  function  z  «  having  A(z)  -  2R  .  Let 


Now  if 


*  »  Yd  -  |n|)«  y  >  o. 
Y  is  chossn  to  achiavs 
with  C2  depending  upon 


Pro*  (2.22)  it  follows  that  Y2^*  <  Ms)  <  Y2^*0*-1** 

A  2 

A( s)  »  2R  ,  a  simple  computation  using  (2.32)  shows  1  <  , 

P(-1). 
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ra» 


3.  ESTIMATES  FOR  THE  EXTENDED  PROBLEM 


In  this  section  some  additional  notation  will  be  useful.  When  x2  <  0  the  expression 
for  a  in  (2.21)  reduces  to  a  function  of  and  p2 


2  2  2  ^  2 

P1  +  P2  P1  P2 

*,pi'p2)  *  5r  Terror + 11  -  V  — r~ 


A  functional  sieilar  to  (2.25),  but  associated  with  8^,  is  defined  by  the  expression 


A(w)  «  /  p^(x2)a(Pw)  . 

°k 


(3.2) 


will  be  the  typical  interior  type,  relative  to  the  variable  x,,.  The  constants  occurring 
in  the  estimates  will  be  denoted  by  C,  possibly  with  a  subscript  or  superscript,  or  in 
the  case  of  a  Holder  exponent,  by  the  letter  ou  These  numbers  will  depend  upon  the 
maximum  density  p(-1),  upon  the  positions  of  the  discontinuities  in  p  in  (1.29),  and 
upon  the  sire  of  p'  where  it  is  continuous,  but  will  be  independent  of  e, 6,  and  the 
period  2k.  The  estimates  also  depend  on  a(x2,p1>p2)  and  its  derivatives  with  respect 
to  p,  and  p2,  but  in  an  inessential  way  (cf.  (1),  section  3).  By  lemma  2.4  we  can  also 
absorb  the  dependence  on  1  into  the  constants  referred  to.  Having  indicated  that  the 


parameters  C  and  a  depend  on  the  given  density  p  we  will  usually  not  display  the 
dependence.  In  general  we  still  suppress  the  parameters  e,  6  and  k.  Throughout  this 
section  6  >  0  and  k  <  •»  at  the  end  we  display  the  e  dependence  and  let  e  approach 
sero.  The  issMdlats  aim  is  to  obtain  estimates  of  w  and  its  derivatives  in  terms  of 


(3.5) 


Integral*  of  )7w|2  over  subregion*  of  8^.  It  follow*  fro*  (2.22)  that 

/  |7w|2  <  C  Mw) 

°k 

and  thus  the  various  norms  of  w  can  be  estiaated  in  terms  of  the  sis*  of  R^  ”  Mw). 


3.1 ,  There  is  an  a  >  0  such  that  the  solution  w  in  proposition  2.2  satisfies 

Iwl2  <  C  /  |7w|2  (3.6) 

C°(8-)  a* 


Proof.  This  is  immediate  fro*  the  known  results  on  elliptic  equation*  (  [25]  ,  Theorem 
8.29),  for  w  satisfies 

K7 


with  atj  -  P  /  »1j(x2't7w)dt' 


Remark.  The  symbol  a  occurring  in  subsequent  results  and  in  Theorem  1.2  should  be 
understood  to  be  the  smaller  of  the  exponents  occurring  in  lesnas  3.1  and  3.2. 


Lemma  3.2.  Let  w  be  the  function  occurring  in  Theorem  2.2  and  let 
is  an  a  >  0  such  that 


and 


IVvl2  <  C  /  |7w|2 

L  (8")  8’ 


Ivl 


c“(B") 


Then  there 

(3.8) 

(3.9) 


Proof.  This  lemma  combines  lemmas  3.1  and  3.2  of  [1],  and  the  proofs  are  essentially  the 
same.  First  one  shows 


/  | 7v|2  <  C'  /  |7w|2  (3.10) 

8"  8* 


for  8"  C  8"  C  8'  using  a  cutoff  function  adapted  to  8"  and  8'  and  a  test  function 


4  ”  C2v  with  (2.27)  (equivalently  (2*29)  with  w+  •  w) .  Inequality  (3.9)  Is  then  a 
3x1 

consequence  of  theorem  8.29  of  (25).  completing  the  proof. 

2 

In  addition  to  a  global  estimate  of  Vv  in  L  the  following  local  estimate  will  be 
used.  For  convenience  we'll  assume  0  <  a  <  V2  in  Lemma  3.2. 


3.3.  Let  v 


wx  be  as  in  the  previous  le 


and  suppose  a(w)  -  R  .  Let 


x  -  (x. ,x„)  be  a  point  in  Q  and  b  C  0  the  ball  of  radius  a  <  1/4  centered  at  x. 
1  2  0 


Then 


(/  |Vv|2)1/2  <  C  R'0a 


(3.11) 


where 


R* 


/  / 


V1 


-1  - 


ivwrdx1dx2 


x,-1 


and  a  is  the  exponent  from  lemma  3.2. 

Proof.  He  use  the  ideas  of  (25).  chapter  12.  Let  a  be  a  radial  coordinate  with  respect 
to  an  origin  at  x  and  let  f  ■  f(o)  >  0  be  1  C  function  with  support  in  B2(j 

(suppose  B2<j  C  Qj  otherwise  extend  w  to  be  odd  and  p  even  about  x2  ■  ±1 ,  obtaining 
a  weak  solution  of  an  equation  on  a  larger  region).  Suppose  that  5  1  on  Bg  and  that 
|Vf|  <  2/o.  Let  h(x2)  «  -IplXj)  and  note  that  |h|  has  a  bound  depending  on  the 
maximum  density  p( — 1 ) .  If  4  is  any  test  function,  (2.27)  yields 

1  PW*,>  dr  ♦  ■  /  (db  (hw)  - h  |r)*  • 

i  2  2 

_  3  2 

With  y  ■  v(x)  and  ♦  »  t —  (41  (v  -  f ) J  the  last  equation  can,  after  integration  by 

3xl 

parta,  be  written  as 

-  /  ^  P*l)  dS;  ^  hw)  37;  -  /  hwx2  37; 


2a 


where  integration  is  over  B, 


unless  indicated  otherwise.  Or 


'  !  C“ik  [2^Xi‘v  -  T)  ♦  »2  §~] 


/  hw  [24-*  (v  -  r)  ♦  V2  v~]  -  /  hw  [2**  (v  -  y)  ♦  *2  -*£-] 

1  2  dX2  *2  1  3xi 


Since  p  >  1,  { h  |  <  c* ;  aik  >  v,  11  1  quadratic  form/  and  <  C*, 

10..1 2 


/  riVvr  <  C  /  2|7v|*|7*|  | v  -  T I 


♦  C'  /  2|7w|*|7*||v  -  rl  *  C'  /  *2(|w  |~|  +  I*  |~|) 


X1JX2 


x  3x 
2  1 


2  12 

Using  the  inequality  2ab  <  Ca  +  —  b  to  abaorb  terns  involving  7v  into  the  left  member 
of  the  last  inequality  and  to  combine  terms  involving  7w  we  arrive  at 


/  f  1 7v |  <  C,  /  |7*|2(v  -  y)2  ♦  C2  /  *2| 


2,7w|2 


(3.12) 


Now  let  C  *  c(x. )  be  a  cutoff  function  which  is  equal  to  1  on  B  and  vanishes  for 
1  2o 

Xj  ■  Xj  -  1 .  Then 


,‘W  "  1  T.  «<“)V(,'x2)d* 


V1 


x1+2o  x2+2o  x1 

II  ♦2*'x  <  /  /  *2(/  (C*w  +  O  )de)2 

B  2  -  —  .  S  S*. 

20  x1-x1-20  x2-x2-2o  x^-1 


x1+2o  x2+20  xf  +  1 

<  C'l  [/  /  <w2  +  c2w2  )  dsdx  ]dx  <  C”(R')2a  (3.13) 

—  —  ■  *  • 
x^-20  x2~20  x^-1  14 

where  use  has  been  made  of  (3.8).  Since,  by  (3.9),  v  »  w  has  a  0°  norm  bounded  by  a 

X1 

multiple  of  R* ,  it  follows  from  (3.12)  that 

/  *2|7v|2  <  C1(|)2(R')2o2°.o2  *  C2(R’  )2(o  +  o2)  . 

It  is  assumed  that  0  <  a  <  V2  so  (3.11)  follows. 

1  N 

Recall  the  notation  n  <  . . .  <  n  for  the  points  of  discontinuity  of  the  density  p 


given  in  (1.29). 
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tj—ji  3.4.  Let  w  be  the  solution  occurring  in  proposition  2.2  and  let  r  be  the  cutoff 
parameter  associated  with  a  in  (2.21).  Then  there  is  a  constant  C  such  that  for 
-1  *  *2  <  0 

Wx  (X1'X2)  <  2rP(T'N)/P(0)  +  CR’  (3.14) 


while  for  -1<x2<-e<0 

"Wx  (x1»x2)  *  2'‘P<n")/p(0)  +  CR’  (3.15) 

where 


R  •  -  /  /  |Vw(s,t)|2dsdt  (3.16) 

-1  xr2 

Proof.  If  -1  <  x2  <  nN/2  the  result  follows  by  using  a  comparison  argument  in  the  region 

O'  -  { (s (t)  |  Is  -  xj  <  1,  x2  <  t  <  0} 

The  fact  that  p*(t)  <  p’(t)  allows  the  proof  of  lemma  3.3  of  [1]  to  be  carried  over  with 
constants  depending  on  the  width  In  I . 

Por  V2  nN  <  x2  <  0,  the  inequality  (3.14)  is  obtained  from  a  comparison  argument  on  a 

N  N 

region  where  the  second  variable  is  between  h  and  Xj.  Let  d  “  |n  -  *2 1 .  To  simplify 
notation  let  (x^Xj)  be  the  origin  of  new  coordinates  (x{,x£)  and  then  omit  the  primes 
to  obtain  a  region 

5  -  {(x1#x2)  |  IXjl  <  1,  -d  <  x2  <  0} 

We  continue  to  use  the  expression  p(x2)  for  the  density  in  the  new  coordinates  and  let 

-q  (q  >  0)  be  a  lower  bound  for  p  or  0.  Note  that  a  in  (2.21)  is  independent  of 

*2 

spatial  coordinates  in  the  region  under  consideration.  Let  Q  be  defined  by 

"  3^  p(x2)\<7w) 

and  observe  that  for  the  solution  w,  Qw  »  Xp’w  <0.  If  u  is  constructed  so  that 

u(0,0)  •  w(0,0),  u  4  w  on  30,  and  Qu  >  0  in  0,  then  according  to  theorem  9.2  of 

[25],  u  4  w  in  0.  It  will  then  follow  that  w  (0,0)  4  u  (0,0).  Por  u  take 

x.  x_ 


*  ~  • 
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(3.17) 


u(x.,x_)  *  w(0,0)  +  w  (0,0 )x,  A  g  ♦  A,e 

1  c  1  Q  J 


where 


(a  fra 


1  +  CI 

g(x1,x2)  -  Re([-(x2  +  iXj ))  ) 


3.2)  and 


e(x_)  »  /  (ad  +  2p(-d)/p(a)  )da 


If  |u  |  is  larger  than  the  cutoff  r  then  (cf.  (2.21)) 
*2 


3  ,  ,  3u 

'  3^  ptx2>  3^ 


p  A  g  ♦  A_d 
x2  a^Xj  3 


-1 


since  g  is  harmonic  in  ft.  Then  since  |g  |<(1+a)(1+d), 

X2 


A3  *  1  +  <*)(d  +  d  ) 

will  insure  that  Qu  >  0  on  ft.  Since 

Ux  *  Aa9x  *  A3(X2d-1  *  2p(-d)  )/p(x2 )  >  -Ao(1  ♦  a)(1  ♦  d2)  +  A3  , 


(3.18) 


having 


A  ~A  (1  ♦  o){1  +  d  )  >  r 
3  a 


(3.19) 


suffices  to  give  |u  I  >  r.  The  term  w(0,0)  +  w  (0,0 )x  is  of  order  R'  according  to 
*2  X1  ,+a 
lemmas  3.1  and  3.2  and  since  g  is  negative  and  of  order  | x ^ |  when  «j  *  0  it  will 


suffice,  as  in  [1],  lemma  3.3,  to  choose  A  >  C'R'  to  have  u  <  w  on  that  part  of  3ft 

01 


where  x2 


As  in  the  lemma  cited,  u  «  w  will  also  be  satisfied  where  Xj  *  11  and 


x.  <  x.  <  0  for  some  x_  e  l-d,0),  with  such  a  choice  of  A  .  On  the  remainder  of  the 
2  2  2  Gl 

mm  ^ 

boundary  of  ft,  i.e.  where  x_  <  x_,  e  <  /  1*ds  <  x.  <  0  Since  the  first  three  terms 

*  1  0  2 

on  the  right  side  of  (3.17)  are  of  order  R',  a  choice  of  A^  >  C*R'  will  make  u  <  0  on 

this  remaining  portion  (recall  w  >  0,  so  u  <  w).  The  coefficient  Aj  can  be  increased 

if  necessary  to  satisfy  (3.18)  and  (3.19),  and  can  ultimately  be  chosen  to  be  of  size 

r  ♦  CR’.  Since  u  (0,0)  ”  2A  p(-d)p  "* ( 0 )  (in  the  new  coordinates)  a  bound  for 
x2 

w  (0,0)  follows)  namely,  the  upper  bound  for  w  contained  in  (3.14).  The  bound  from 
x2  X2 
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below  follows  in  a  like  Banner,  using  the  fact  that  pj  £  “  p'^  for  <  -e,  and  hence 
can  be  estimated  independently  of  e. 

Le— a  3.5.  Define  6  by 


0( e,v,w) 


/  (iT  6  >  0 

Qk 

/  v(x1 ,0)w(x,0)dx1 ,  e-0 


(3.20) 


where  t£  is  the  function  occurring  in  (2.19)  and  »  {(x^.O)  |  |x^ |  <  k}.  Suppose 
•*  0,  w  in  H^,  and  -*•  v  in  for  J  »  1,2,...,  the  convergence  in 

being  weak.  Then 

lia  0(c  ,w  ,v  )  -  8(0,w,v)  (3.21) 

jam  3  3  3 

Proof.  An  integration  by  parts  produces  the  expression 

8<c,v,w)  -  8(0, v,w)  -  /  T£(v  w  *  vw  )  .  (3.22) 

°k  2  2 

1/2 

There  in  n  continuous  linear  map  taking  an  element  of  to  its  trace  in  H  (1^) 

2 

((26],  theorem  9.4).  Since  this  last  space  embeds  compactly  into  L  ( ) ,  the  trace  of 

2 

Wj  (or  Vj)  converges  strongly  in  L  (T^).  The  weakly  convergent  elements  w^  and  v^ 

3  3 

lie  in  a  bounded  set  in  Hv,  from  which  it  follows  that  f-r-—  v  V  +  v.  -r— -  w.  lies  in  a 

K  J  J  J  ^ 

bounded  set  in  Lp(fl  )  for  any  p  <  2.  Since  T  converges  to  zero  in  Lr  for  any 
k  ej 

r  <  «•,  it  follows  from  Holder's  inequality  that  the  integral  term  in  (3.22),  evaluated 
with  e  -  c^,  v  •  Vj,  and  w  •  w^,  converges  to  zero  as  j  ♦  ». 

Proposition  3.6.  For  each  S  >  0,  k  >  0  and  R  >  0  there  is  a  X  >  0  and 

w  «  h£(Q)  n  C2(0>  n  C0,1(0)  with  w  >  0  and  a(w)  ■  R2  such  that 

a’(w)  -  Xb' (w) 

in  the  notation  (2.18)  and  (3.2).  That  is 
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/  p-a  (Vw)$  “  X  J  plw$  (3.23) 

°x  Xi  °k 

for  all  4  €  H^(fl).  The  function  w  is  an  extremal  for  the  problem 

i  *\2 

sup  (-1)  /  Pl(x  >  (3.24) 

A(u)-R  St 

u«H*(Q) 

Moreover,  with  the  restrictions  that  B  C  fl  and  that  fl  is  replaced  by  fl  in 

0 

definitions  (3.3)  and  (3.4),  the  estimates  (2.30),  (3.6),  (3.8),  (3.9),  and  (3.11)  hold 

uniformly  in  6  and  k,  as  well  as 

Iw  I  <  2rp(hN)/p(0)  ♦  CR’  (3.25) 

X2  L“(fl> 

with  R'  defined  as  in  (3.16). 

Proof.  Let  wE  e  H*(Q)  denote  a  solution  obtained  from  proposition  2.2  for  0  <  e  < 

and  let  flE  •  {(x^.x^)  e  fl  I  <  -e}.  From  (3.5)  and  lemma  3.1  we  know  that, 

independently  of  6,k  and  e,  the  functions  w£  lie  in  bounded  sets  in  H^(fl)  and 

C°(fl).  Likewise  from  lemmas  3.2  and  3.4  the  derivatives  wE  lie  in  a  bounded  set  in 

X1 

a  A  £  M  £ 

C  (Q)  and  the  derivatives  w  are  bounded  in  L  (Q  ).  In  addition,  since  p,  is  in 

x,  6,  t 

1  fl  ^ 

C  *  (-1,0),  it  follows  from  quasilinear  elliptic  theory  U24),  theorem  6.3,  p.  283)  that 
for  each  n  >  0,  wE  e  C2'^(flr')  with  bounds  depending  on  i,  but  uniform  in  E  as  c 

approaches  0 . 

Let  e  take  the  values  2”1,  j  -  1,2,...,  and  let  w^  denote  the  corresponding 

solution.  By  the  Arzela-Ascoli  theorem  a  subsequence  w1 ^ ,w^  2, • •  •  converges  in  c°(6) 

2  1/2 

and  in  C  (fl  )  to  a  function  w.  A  further  subsequence  wj  ,, w2  2,...  converges  in 

2  1/4 

C  (fl  )  to  w.  Continuing,  we  find  the  usual  diagonal  sequence  wj  “  wj  j  converging 

to  w  in  C2(fln)  for  each  n  >  0.  Moreover,  the  estimates  cited  at  the  outset  of  the 

proof  give  bounds  on  w  8  C°(fl)  and  on  -1^-  e  C°(fl)  as  well  as  the  estimate 

3x1 

Iw  I  m  <  2rp(nN)/p(0)  +  CR’  . 

2  L  (fl) 
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2  *  o  1  • 

It  follow*  that  wee  (fl)  and  ha*  an  axtanaion  to  fl  which  la  In  C  '  (fl).  Clear ly, 
than,  w  e  H*(Q). 

Next,  It  follow*  from  13.5)  that  Iw .  I2  <  CA(w  )  -  CR2  and  «o,  for  a  further 

3  "k  3 

subsequence  (alao  denoted  w^),  it  can  be  assumed  that  w^  converge*  to  w 

weakly  In  (fl)  and  atrongly  In  Lj ( fl^ ) .  Sine*  a(p^ ,p2 )  la  a  convex  function,  th* 

2 

functional  Mw)  i*  alao.  Then  the  aet  of  u  for  which  Mu)  <  R  ia  convex  and,  aa 

w3  converge*  weakly,  A(w)  <  R2.  Let  p,  with  c  ”  e .  be  denoted  by  p.,  let 

o,  e  3  3 

2b^  “  /  Pj(Wj)2,  and  let  2b  ■  /  p’w2  (here  and  in  the  remainder  of  th*  proof  all 
integral*  ara  taken  over  0^  unless  otherwise  indicated).  It  ia  then  immediate  from  the 
known  convergence  of  w^  and  lemma  3.5  that  b^  ♦  b  a*  J  ♦ 

Let  the  extreme  value  in  (3.24)  be  denoted  by  b.  First  we  ahow  that  b  <  b  is 
impossible,  if  b  -  b  “  d  >  0  then  from  the  characterization  (3.24)  there  ia  a 
u  e  H*(fl)  satisfying  A(u)  -  R2  and  /  Pj(u+)2/2  >  b  -  d/3.  Then  from  lemma  3.5,  for 
all  sufficiently  small  positive  e,  /  pi  (u+)2/2  >  b  -  d/2.  But  then  since  wE  ia  an 

o,  E 

extremal  for  th*  problem  with  £  >  0,  /  p^  £(we)2/2  >  b  -  d/2.  This  is  incompatible  with 

m.  2 

having  bj  converge  to  b  “  b  -  d  as  j  ♦  «*.  We  have  shown  that  A(w)  <  R  and 

/  Pjw2/2  >  b.  However,  neither  of  these  inequalities  can  be  strict  without  contradicting 

the  characterization  (3.24).  In  particular,  if  A(w)  <  R2,  one  easily  shows  that 
2 

A(tw)  ■  R  for  aome  t  >  t  with  a  corresponding  supremum  larger  than  b  in  (3.24). 

e 

As  regards  the  equation  (3.23),  since  w  is  even  in  x2,  we  conclude  from  (2.27) 
that  for  each  e  >  0 

/  Pj(x2 )*1(7we)  |£-  -  A 6  /  p^£w%  (3.26) 


for  all  +  e  H^tfl).  It  is  easy  to  verify  that  for  i  •  1  or  2, 

la1(p1,p2)|  <  C(|p,l  ♦  |p2l) 

(cf.  [1]  lemma  2.1)  and  thus  for  4  *  C0,1(flk)  the  integral  of 
will,  by  Holder's  inequality,  be  Otn^2)  uniformly  as  c  ♦  0. 
left  hand  member  of  (3.26)  converges  to  the  left  hand  member  of 


(3.27) 

p«*i\  OV,r  °k  * 

Thus  with  t  m  Cy  ths 

(3.23)  as  J  ♦  «*.  Since 
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X  ^  satisfies  the  inequalities  (2.30)  it  can  be  assumed,  without  loss  of  generality,  that 

X  J  converges  to  a  value  X  satisfying  (2.30).  it  then  follows  froa  lenu  3.5  that  the 

right  hand  member  of  (3.26)  converges  to  the  corresponding  term  in  (3.23).  Thus  (3.23) 

holds  for  4  «  C0'1  and  the  extension  to  (6  1^  follows  by  continuity. 

He  next  show  that  (3.8)  holds  for  the  w  just  obtained.  Let  v  ”  (w  )  and 

J  J 

v  •  .w  e  From  (3 #8)  it  follows  that  for  y  >  n  >  0 

^e 


/ 

Q-nft 


Y 


'1  '  n  |VV 
1  Q*non  3 


♦  c,  /  n  •’-j' 

o'\an  3 


and  hence  for  each  Y  >  0, 

/  y  |Vv|2  <  C.  /  |Vw|2  (3.28) 

OTOt  ft1 

provided 

_  It  0 

lia  lia  /  J  |7w  |2  -  0  (3.29) 

n+o  j-Mt  -V  -n  3 


If  (3.29)  does  not  hold,  then  froa  (2.22)  and  the  observed  convergence  of  A£  to  A, 
uniformly  on  bounded  seta  in  Hk,  it  will  follow  that  there  is  a  sequence  ♦  0, 
a  -  1,2,...  such  that 

_  k  0 

lia  /  /  pa(*r  )  >  (X V  >  0  (3.30) 

j*-  -k  -nm  3 

for  all  a.  But  then 


k  a  k  a  2+2 

J  J  pa(Vw)  -  lia  /  /  pa( Vw . )  <  R  -  (R  )  (3.31) 

-k  -1  -k  -1  3 


2 

for  all  a.  If  a  approaches  infinity  (3.31)  yields  A(w)  <  R  ,  which  has  been  ruled 

out.  It  follows  that  (3.28)  holds  for  each  Y  >  0  and  hence  with  the  integral  taken  over 

all  of  fl",  completing  the  proof  of  (3.8).  The  other  inequalities  listed  in  the 

proposition  are  done  similarly. 

2 

Since  w  e  C  (fl)  the  following  result  is  iaaediate  froa  (3.23)  and  the  strong  aaximua 
principle  (cf.  (3.7)). 
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4.  A  RETURN  TO  THE  ORIGINAL  PROBLEM 


Here  we  show  that  the  solution  (A,w)  of  A*(w)  -  Ab*  (w)  given  by  proposition  3.6 
is,  for  suitably  restricted  R,  a  solution  of  the  "physical*  problem  (2.12).  We  also 
obtain  additional  estimates  which  complete  the  assertions  of  theorem  2.2  (except 
exponential  decay)  in  the  case  4  >  0  and  k  <  “. 

In  section  3  estimates  were  derived  for  w  and  v  ■  w  .  The  norms  of  w  and  v  in 

a  X1 

C°  and  of  Vv  in  L2  on  a  region  fl"  C  Q  were  estimated  in  terms  of  R',  the  L2 

norm  of  Vw  on  a  larger  region  fl*  C  fl  (cf.  proposition  3.6).  In  addition  in  (3.11) 

the  L2  norm  of  Vv  on  a  ball  B  C  fl  was  shown  to  be  of  order  Roa  and  in  (3.25)  an 

O 

L  bound  was  given  for  w  in  terms  of  r  and  R*.  Our  next  step  is  to  obtain  bounds 

x2 

on  w  in  terms  of  R*  alone.  We  require  a  preliminary  lemma.  Recall  that  A  in  (3.2) 
2 

involves  a  from  (2.21)  with  a  cutoff  parameter  r. 


Lesssa  4.1.  Let  w  be  the  function  occurring  in  Proposition  3.6.  There  exist  positive 

constants  r^  and  Rg  such  that  if  A  is  defined  using  a  cutoff  r  <  r^  and 
2  2 

A(w)  «  R  <  R  ,  then  v  -  w  satisfies 

U  X, 


I  l?v  |2  <  C  /  |Vw|2 

S"  xi  s* 

where 

S"  -  {(x^Xj)  |  b(  <  <  b2,  -1  <  x2  <  nN/2}  , 

S'  “  {(x1#x2)  |  bj  -  1  <  Xj  <  b2  +  1,  -1  <  x2  <  nN/4}  , 

s*  c 


(4.1) 


(4.2) 


Proof.  This  lemma  is  a  counterpart  to  lemma  3.4  of  [1]  and  the  proof  is  similar.  One 

a  2 

starts  with  a  test  function  4  -  - —  5  v  in  (3.22).  Here,  however,  the  cutoff  function 

3x1  X1 

is  taken  to  be  1  on  S"  and  to  vanish  outside  of  S*  given  by  (4.2).  Now  that  5 
depends  on  x2  as  well  as  on  x1 , 


2  _3 

ax 


x  *  2«x  wx 
2  1  2  1 


In  th«  earlier  caae  just  cited,  the  term  with  t  was  not  present,  but  the  new  term  is 

X2 

N 

easily  incorporated  into  the  estimates  in  terms  of  a  constant  involving  the  gap  Ini. 


Lemma  4.2.  Let 
constant  r  < 

where 


w  and  Rg  be  as  in  the  previous  lemma.  Then  there  is  a  positive 

such  that  if  A  is  defined  using  r  «  r,  and  0“  C  Q‘  C  (2^, 

Iw  I  <  CR*  (4.3) 

X2  L  <Q*> 


R’  -  (/  |Vw|2)1/2 

O' 


(4.4) 


Moreover,  if  <2  is  a  closed  region  not  containing  any  points  of  discontinuity  of  p  in 
(1.29)  then 

Iw  I  <  C'R*  (4.5) 

*2  C  (OmS) 


for  all  sufficiently  small  &  where  a  is  the  exponent  from  lemma  3.2. 

U 

Proof.  The  estimate  (4.3)  for  x^  <  n  /2  is  shown  as  in  lemma  3.5  of  [1]  with  one  small 

change.  In  the  present  context  the  function  w  need  not  vanish  on  each  line  x,  - 

*2 

constant.  However,  since  Iwl  <  C'R',  w  is  of  order  R*  at  some  point  (x  ,x  )  for 

L*  X2  2 

each  x1 ,  by  the  mean  value  theorem,  and  this  suffices  for  the  argument.  An  argument 

given  in  the  lenmaa  cited  also  shows  that  the  oscillation  of  w  over  a  distance  d 

*2 

satisfies 


osc(w  )  <  CR'd 
2 


(4.6) 


in  the  region  where  x  <  1/2,  where  C  does  depend  on  a  bound  for  |p  I  in  the  region 

*2 

(and  hence  is  independent  of  6,  for  small  4,  in  a  region  $2  containing  no 
discontinuity  of  p).  In  the  present  context  lemma  4.1  is  used. 

How  consider  the  region 

Q°  -  {(x1,x2)|nN/2  <  x2  <  0} 
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and  recall  that  IP  I  ia  bounded  In  Q  ,  Independently  of  6.  Proa  Corollary  3.7,  w 
2 

satisfies  the  equation  (3.31)  In  £T  i  that  la, 

P(Xj)  I  a^lVw) w^  ^  +  p'SjfTw)  -  Xp’w  .  (4.7) 

0 

If  C  (1  n  Q"  l*  any  disc  of  radius  a,  then  frca  (4.7)  it  follows  that 


/  lw  I2 

4  1  V  * 


/  “T  ^  •ijwx  X  "  p”  42(W>  *  X  "l*  (4*8 

B_  af_  2<i+j<3  3  11  p  *  p 


•  ■  v  jp  •  2  t  «• 

Bo  22  Bo  *22  2<i+j<3 


The  expressions  la^/Sjjl  an4  I P* /P I  ®r®  bounded  above  in  Q°  so  the  expression  on  the 
right  in  (4.8)  can,  according  to  proposition  3.6,  be  estimated  using  lemmas  3.1  to  3.3 
together  with  inequalities  (3.13)  and  (3.27).  The  result  is 


/  |w  I2  <  C(R')2(02“+  o)  (4.9) 

B  2  2 
0 

where  a  (assumed  to  satisfy  0  <  a  <  V2  )  is  from  lemma  3.2.  The  estimate  (4.9)  together 


with  lemma  3.3  shows  that 


/  |Vw  |2  <  C(R*)V“ 


4  •  x  ■  -  - 

Bc  2 

and  hence  by  Horrey's  estimate  ((25),  p.  158,  combined  with  the  Schwarz  inequality)  the 
oscillation  of  w  over  a  distance  d  in  0°  satisfies  (4.6).  This  estimate  in  fl° 
combined  with  the  estimates  for  <  n  /2  yield  the  assertions  of  the  lemma. 

Prom  lemmas  3.1,  3.2  and  4.2  the  following  result  is  immediate. 


Corollary  4.3.  Under  the  conditions  of  the  previous  lemma 

IwB  (tl)  <  CR' 

C 0,1 


(4.10) 


Iwl  <  C'R' 

c  '  (0"nfl) 


■**  ? 


-  V 


In  what  follows  the  notation  (2.18)  for  F(w)  and  B(w)  will  be  used  and  for  the 


renal nder  of  this  section  an  Integration  Is  over  SI  unless  otherwise  indicated. 


Proposition  4.4.  There  la  an  R^  >  0  depending  on  the  density  p  In  (1.29)  such  that  for 
each  8  in  (0,8^),  k  >  0,  and  R  In  (O.R^ )  there  Is  a  X  >  0  and  a  nonnegative 
w  e  h£(S))  n  c2(Q)  n  c1,a(5)  with  p(w)  -  r2  and 

P'(w)  -  Xb'(w)  .  (4.12) 

The  pair  (X,w)  is  a  solution  of  the  problem  (2.12),  i.e. 

plXjJf^Vw)  -  Xp'w  In  S^  (4.13) 

and 


w  "  0  on  x2  »  -1 

p(Xj)f2(Vw)  ■  Xw  on  *2  “  0 


(4.14) 


where  P  m  Pg  (cf.  (1.34)).  The  estimates  (4.10)  and  (4.11)  from  corollary  4.3  hold 
uniformly  in  6  and  k. 

Proof.  From  lemmas  3,2  and  4.2  there  is  a  positive  R.,  such  that  for  R  <  R, ,  |fe|  <  f. 

Thus  A(w)  *  F(w)  and  a^Vv)  *  f^Vw)  for  1  »  1,2  so  that  (4.12)  follows  from 

Proposition  3.6  and  (4.13)  from  Corollary  3.7.  Naturally  w  «■  0  for  x2  -  -1  and  all 

estimates  for  X  and  w  persist.  In  particular,  from  Corollary  4.3  it  follows  that  VW 

a 

has  a  uniquely  defined  C  extension  to  the  line  x2  “  0.  we  assume  the  extension  is  made 
and  can  then  explicitly  express  the  boundary  condition  which  is  implicit  in  the  weak 
equation  (4.12). 

Let  g  *  g(x1 )  be  an  arbitrary  2k  periodic,  C  function  and  for  0  <  s  <  1  let 
H(Xj)  -  max{1  ♦  s  ^J,!)}.  Using  (4.12)  to  write  <F*(w),8>  -  X<B '(*),♦>  with 
♦  “  g(x1 JHtXj)  one  obtains 

k  0 

/  /  P(x  )(f  (Vw)g  H  ♦  f  (7w)gs_1) 

-k  -s  1 

k  0 

“  X  /  [  /  p'  (x2 JwgHdXj  ♦  vfXj.OlglXj)  ]dx?  . 


Sine*  p(x2>  approach**  1  aa  x2  approach*#  0  and  p'  1*  bounded  naar  x^  *  0, 
tha  limit  a*  ■  ♦  0  give* 

k  k 

/  f2(Vw(x1,0))g(x1)dx1  -  X  /  w(xf ,0 ) gtXj )dx  . 

-k  -k 

Sine*  g  is  arbitrary,  the  remainder  of  condition  (4.14)  follow*  and  the  proof  is 
complete. 

Tha  material  in  • action  4  of  (11  *how*  that  in  the  case  of  a  fixed  upper  surface,  the 
streamline  displacement  w  can  be  assumed  to  be  even  in  x  for  each  h  and  nonincreasing 
in  x  as  x  runs  from  0  to  k.  That  is,  an  extremal  function  w  tor  the  appropriate 
variational  principle  can  be  replaced  by  its  symmetrisation  (decreasing  rearrangment)  w 
without  destroying  the  extremal  character  (cf.  (1),  (27]).  Here  again  the  proofs  carry 
over  with  soon  minor  modifications  which  are  furnished  by  the  following  result.  The 
functional  N  defined  below  is  part  of  the  functional  4  in  (2.1)  and  enters  in  the  proof 
of  lensu  4.3  of  [1].  Using  evenness,  it  will  suffice  to  consider  4  <  0. 


Lemma  4.5.  Let  u  »  u(x,4)  be  a  nonnegative  piecewise  linear  function  defined  on 

«  l-k,k]  x  ,0]  and  suppose  u  has  a  continuous  extension  to  R  *  [  ,0 ]  which 
is  2k  periodic  in  x.  Let  y(x,4)  -  $(♦)  +  u(x,4)  where  y(  ♦)  is  the  inverse  funtion 
to  (1.16)  and  assume  y^  >  0  a.#..  Define 


N(y) 


1  +  *3 


dxd* 


(4.15) 


so  that 


<H'(y),u>  -  /  dxdf  (4.16) 

°k 

Then  if  u(x,4)  and  y(x,t>  are  the  symmetrisation*  of  u(x,f>)  and  y(x,4), 
respectively,  it  follows  that 

<*!'($), u>  -  <H'(y),u>  (4.17) 


-36' 


and 


(4.22) 


/  (1  +  \  t2|V*|2  +  0(t2)l  4  /  [1  ♦  \  t2 1  V<>| 2  +  0(t2))  . 

s.  s 

?  y 

Raplacing  u  by  y  in  (4.19)  one  sees  that  Sy  and  S^  have  equal  areas.  Thus  the 

contribution  of  the  "1*  in  each  integrand  of  (4.22)  can  be  omitted.  If  the  remaining 

1  2 

inequality  is  divided  by  —  t  and  t  is  allowed  to  approach  zero,  the  inequality  (4.21  ) 
is  the  result. 


Proposition  4.6.  There  is  a  function  w  satisfying  the  conclusions  of  proposition  4.4 
which,  in  addition,  satisfies  w  ”  0  and  w  >  0  where  -1  <  <  0. 

Proof.  As  noted,  the  assertion  that  one  can  take  w  ■  9  in  proposition  2.2  and  thereafter 

follows  from  the  arguments  of  [1],  section  4  together  with  lemma  4.S.  The  positivity  of 

w  in  follows  from  corollary  3.7.  On  we  have  w  nonincreasing  for  0  <  x  <  k 

and  f,(Vw)  «  Xw.  Since  w  (0,k)  "  0,  having  w(0,k)  "  0  would  entail  f  (0,w  )  »  0 

2  X1  2  *2 
at  (0,k).  Since  f2  “  a2  for  the  solution  w  it  follows  from  (2.22)  that 

w  (0,k)  “  0,  violating  the  strong  maximum  principle. 

*2 

The  linearization  of  (4.13),  (4.14)  about  w  ■  0  is  the  problem 

3^  p(V  !^“  Xp,(* 2)z  in  «k 


z  •  0  on  x^  ”  “1  V 


(4.23) 


on  x2  -  0 


The  lowest  eigenvalue  (denoted  |i)  for  (4.23)  can  easily  be  shown,  through  separation  of 
variables,  to  be  positive  and  to  correspond  to  a  function  of  Xj  alone,  £  »  £(x2>  >  0 
which  solves  the  Sturm-Liouville  problem  obtained  by  omitting  x1  dependence  in  (4.23). 
It  also  provides  an  extremal  for 


sup 


1  7  7 

[/  p'r  +  r<o>]  , 

0 


(4.24) 
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the  analogue  of  0.24).  The  solution  i  can  be  obtained  by  a  shooting  method  or  by  using 
(4.24). 


L— a  4.7.  The  "lowest*  eigenfunction  £  for  (4.23)  satisfies  C'lXj)  >  0. 

Proof.  One  integration  of  the  ordinary  differential  equation  for  £  -  £(Xj>  corresponding 
to  (4.23)  gives 

P(O)CIO)  -  plXjJC'lXj)  <  o  . 

since  5* (0)  -  wt(0)  >  0,  C’(x,)  >  0. 

0 

Since  lemma  4.7  implies  /  p£J  >0.  it  is  possible  to  use  a  trial  function 
-Cjxl  -1  *2 

s  »  C0t(x2)e  in  (3.24).  as  in  (1).  section  S,  to  show  the  following. 


Proposition  4.8.  There  are  positive  constants  R  and  k, (R)  such  that  for 
0  <  6  <  6q,  0  <  R  <  R  and  k  >  k1(R)<  the  pair  (X,w)  in  proposition  4.3,  chosen  in 

conformity  with  proposition  4.6,  satisfies 

X  <  U(1  -  CR4/3)  (4.25) 


and 


Iwl  >  C-R4/3 

m 

L 


(4.26) 


where  C  and  C*  depend  upon  p  in  (1.29) 

Note  that  the  inequality 

/  |Pw|2  <  C  max  w2  (4.27) 

0"  O' 

which  is  used  in  obtaining  (4.26)  from  (4.25)  is  valid  in  the  present  context  and  shorn 
exactly  as  in  lemma  3.6  of  (1]. 


Remark  4.9.  The  results  up  to  this  point,  in  particular  propositions  2.1,  4.3,  4.4  and 
4.6,  contain  the  assertions  of  theorem  1.2  (except  exponential  decay)  for  the  case  in  which 
density  "transitions"  take  place  over  intervals  of  width  4  >  0  and  in  which  there  is  a 
finite  period  2k  in  the  horisontal  direction,  all  estimates  obtained  are  uniform  in  5 


and  k  for  aach  "energy"  R  >  0  and  ao  tha  Halting  casa  of  discontinuous  danslty 
{4  ■  0),  Including  tha  praaaura  condition  In  definition  1.1  (11),  follows  by  arguments 
strictly  paralleling  those  given  for  theoraa  8.2  of  [1]. 


S.  EXPONENTIAL  DECAY  AND  SOLITARY  HAVES 

The  ala  in  thia  aaction  la  to  show  that  w  obtained  from  propoeitiona  4.4  and  4.6  haa 

exponential  decay  on  0  <  x  <  k,  aa  doea  the  gradient  of  w.  The  estimate  will  involve 

L**  bounds  on  the  gradient  and  we  reaark  that  in  order  to  get  aatiafactory  estiaatea  for 

w  up  to  x2  •  0  we  had  first  to  pass  to  a  limit,  letting  e  approach  zero  in  the 
*2 

O  0 

expression  £,  here  denoted  p£.  While  there  is  more  than  one  way  to  exhibit 
exponential  decay  we  find  it  convenient  to  reintroduce  the  extended  domain  ft,  the 
function  p£  from  (2.19)  and  the  expression  a  from  (2.21). 


5.1.  Let  v£  be  the  lowest  eigenvalue  of 


3^  p  “  rpcZ'  "  «  x2  «  1 


z(-1)  -  z(1)  -  0 


(5.1) 


then 


We  >  V  -  Ce 


1/2 


(5.2) 


where  p  is  the  lowest  eigenvalue  of  (4.23). 

Proof.  Express  the  odd  function  p£  as  p  ♦  T£  on  -1  <  x2  <  0  as  in 

(2.19)  and  let  €(x_)  >0  be  an  eigenfunction  for  the  problem  (4.23)  corresponding 

1  , 

to  p,  normalized  so  that  /  p ( x2 ) ( C ' )  “  1.  Then  since 
p  >  1,  /  U1)2  <  1,  1 5(  x2 )  -  5(0)|  <  (/*2  (5’)2)V2|x2|1/2,  and 

“  4(0)1  <  .  Extend  5  to  be  even  in  Xj.  The  variational 

characterization  of  p£  enables  us  to  conclude  that 


f  >  -  [  /  p'52  ♦  /  *;e2] 

-1  -1 

>  -  [  /  P’C2  ♦  /  t^(52(0)  -  C,e1/2)]  -  p  -  CC1/2 


-1 


proving  (5.2) 
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If  X  la  fixed  with  X  <  p,  then  for  all  sufficiently  small  e  >  0,  X  < 
according  to  (5.2).  Estimates  involving  the  eigenfunctions  wn>  n  “  1,2,...,  and 
corresponding  eigenvalues  Yn  for 


d  dc  ,1, 

~  P  ~  >  x?^  -  tp« 

2  2 


s(1)  *  *(-1)  -  0 


<5. 3) 


are  done  exactly  as  in  lemma  7.1  of  [1],  merely  by  using  the  fact  that  t P£l  <  P  for 


-1  <  Xj  <  1.  The  Green's  function  for 


L 


t 

X 


Xp^a 


(5.4) 


with  zero  boundary  conditions  at  x2  •  ±1  has  the  form 
g(x  -  x*,n,n')  -  l  * 


y’/S|x~x'| 


n-1  27 


V2 - V^V'1 


(5.5) 


with  7,  >  c0<)*€  -  X). 


Proposition  5.2.  The  function  w  from  proposition  4.6  satisfies 


|w(x,n) I  <  Ca 


8x 


(5.6) 


|7w(x,h) 1  <  C's 

1/2 

for  0  <  x  <  k,  for  any  $  satisfying  0  <  $  <  Y^  .  The  constants  C  and  C*  depend 
on  P  in  (1.29),  R,  and  8. 

Proof.  Define  V  «  (V^.Vj)  where 

Vil*2,Pl'P2)  “  P<pi  ”  ■i<x2,p1'p2>)  <5,7) 

for  i  «  1,2  and  aA  *  3a/3p^  as  befora.  The  function  has  the  mama  parity  as  aA 

(cf.  the  proof  of  proposition  2.2);  i.e.,  a,  is  odd  in  p1  and  even  in  (Xj.Pj), 

2  2 

while  Vj  has  the  opposite  parity.  Both  v1  and  V2  are  of  order  +  p2  near 

2  2  ~2 

P1  *  p2  •  0  and  vanish  for  p(  ♦  )  2r  , 
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Suppose  w  fro*  proposition  4.6  ia  extended  to  ba  avan  on  0  and  than  axtandad  to  ba 
odd  about  Xj  “  ±1.  Convolution  of  w  with  raapact  to  a  mollifying  karnal  which  ia 
symmetric  in  x1  and  x2  and  haa  support  in  tha  ball  of  radiua  t  at  (0.0)  will  produca 
a  family  of  c  functions  *£  which  ara  avan  in  x2,  aatisfy  «(  •  s  ,  and  convarga 
to  w  in  aa  c  approachaa  zaro.  Sinca  |Vw|  <  r,  it  can  alao  ba  assumed  that 

7.3  of  Ml  it  ia  aaaily  assn  that 


|Vs(|  <  r.  Aa  in  Is 


CR 


2/3 

l*e(x,,x  >|  < 

X1 


(5.8) 


on  0  <  Xj  <  k  whara  C  dapanda  on  p  from  (1.29). 


Now  lat  w  e 


ba  tha  waak  solution  of 
c  e 


LXW 


div  V(x2,7«c)  . 


(5.9) 


Given  tha  paritiaa  of  V1  and  V2  it  ia  aaay  to  varlfy  that  w^(x1 ,x2 )  -  we(-Xj,x_ )  is 

also  a  waak  solution  of  (5.9)  and  sinca  LE  ia  coarciva,  w *  ■  w C>  i.a.  wE  is  avan  in 
c 

Xy  Likawisa  w  ia  avan  in  x2  and  2k  periodic  in  x,.  Por  j  “  0,1,2,. ..,2k  -  1 
dafina 

b,  -  /  |7w*I2»  b  -  /  |Vxl2  (5.10) 

3  B  3  B 

J  3 


where  ■  {(x^,x2>  e  a  |  j  <  x(  <  )  ♦  1),  For  a  positive  integer  n  laaa  than  k  let 

b  denote  the  sequence  {b^}  and  b,  tha  saquanca  {b^ }  for  n<j  <  2k  -  n  -  1.  Now 

1/2 

if  6  is  a  raal  number  satisfying  0  «  0  <  ,  it  follows  aa  in  tha  proof  of  lease  7.4 

of  [1],  (one  can  smooth  tha  sign  urn  function  in  a(x2,p^,p2)  and  then  pass  to  a  limit)  that 

bj  «  t ^(b)  (5.11) 

for  n  +  1  <  J  <  k  where 


2k-n-1 


Tj(b)  -  c^b^,  +  bd  +  bj+1)2  +  c2[^  l  .-«1-3lbi  ♦  q,~«3-n>]  <s,12) 


and 
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t  v  ^  X'  — 


22: 


Tha  solution  up  taking  r  (  to  w*  ia  eoftti»w«xi*  from  H^tO)  to  H^(0) 
e  ♦  0  and  i£  convargaa  to  w  in  H^,  w*  converges  to  a  function  w°  in 
weak  fora  of  (5.9)  raquiraa 


and  thua  aa 

V  Th* 


J  [p  1^7 1^7 +  XpV*“ 1  Vv7**1  iJ- 

fik  1  1  Sc 

for  all  ^  e  Hk<fi),  where  pf  la  atill  assumed  to  ba  odd  in  Xj.  Now,  lotting  e  ♦  0 
and  using  laaua  3.5  to  dafina  /  ^w°  ♦  aa  tha  Unit  of  ]  pS<ef,  ona  concludes  that 


But  froa  proposition 
( extended  evenly  in 


1  l»  Ir  ii-  *  ‘  »S«M  ■  1  w->  th 
\  \ 


3.6  it  follows  that  (5.13)  is  satisfied  if  w 


Xj).  Letting  wu  -  w  -  z  it  follows  that 


(5.13) 

is  replaced  by  w 


/  p|V«|2  ♦  Ip's2  -  0  .  (5.14) 

Sc 

Since  the  coerclvity  of  ia  uniform  for  all  small  e,  it  follows  that  s  «  0.  That 

is,  -E  and  s  both  converge  to  w  in  Hb  as  e  approaches  0.  But  then 

b}  <  T^(b)  (5.15) 

where  b^  is  now  associated  with  w°  ■  w. 

The  inequality  (5.15)  will  imply  there  is  an  n,  independent  of  k,  such  that 
b^  <  Ce  2®^  for  n  «  j  <  k  according  to  lesna  7.2  of  [1],  provided  the  sequence  bj 
meets  two  other  conditions.  The  first  is  a  symmetry  condition  b^+i  -  ^k-i-i  (this  shift 

in  index  should  appear  in  (1)>  tha  proof  is  almost  identical)  which  follows  from  tha 

-2/3 

avennese  and  periodicity  of  w.  The  second  ia  a  decay  b^  <  CQ j  which  follows  from 
lemmas  3.6  and  7.3  of  [1],  trivially  adapted  to  tha  present  context.  The  exponential  decay 
of  bj  implies  exponential  decay  of  w  and  Tv  according  to  lemma*  3.1,  3.2  and  4.2. 

This  completes  tha  proof. 
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Remark  5.3.  The  resolution  of  problem  S  0.32),  tha  solitary  wave  problaa,  la  carried 


out  by  letting  tha  period  2k  approach  infinity  just  as  in  theorems  8.3  and  8.3*  of  {1} 
for  3-0  or  i  >  0,  respectively.  The  assertions  of  theorem  1.2  concerning  exponential 
decay  and  the  case  k  »  are  thus  covered,  completing  the  proof  of  the  theorem. 


REFERENCES 


1.  R.  E.  L.  Turner,  Internal  waves  in  fluids  with  rapidly  varying  density,  Annali  della 
Scuola  Normals  -  Pisa,  Ser.  IV,  S  (1981),  513-573. 

2.  J.  L.  Bona,  0.  K.  Bose,  R.  E.  L.  Turner,  Finite  amplitude  ateedywaves  in  stratified 
fluids.  Mathematics  Research  Center  TSR  No.  2401,  University  of  Wisconsin,  Madison, 
1981 i  Jour,  de  Math  Pure  et  Appl.,  to  appear. 

3.  J.  V.  Nehausen  and  E.  V.  Lai tone.  Surface  Haves,  Handbuch  der  Phyaik,  Vol.  IX, 

StrUmungsmechanik  III  (Ed.  S-  FMlgge)  Springer-Verlag,  Berlin,  1960. 

4.  C.  J.  Amick  and  J.  F.  Toland,  On  finite  amplitude  solitary  water  waves.  Arch.  Rat. 
Mach.  Anal.  76  (1980),  9-95. 

5.  A.  M.  Ter-Krikorov,  ThSorie  exacts  dee  ondes  longues  atationnairee  dans  un  liquids 

hdtdrogene  J.  Mdcanlque  2  (1963),  351-376. 

6  K.  0.  Friedrichs  and  D.  H.  Hyers,  The  existence  of  solitary  waves.  Comm.  Pure  Appl. 
Math.  7  (1954),  517-550. 

7.  M.  A.  Lavrentiev,  A  contribution  to  the  theory  of  long  waves,  Aaer.  Math.  So c. 
Translation,  No.  102,  Providence,  RI,  1954. 

8.  P.  R.  Garabedian,  Surface  waves  of  finite  depth.  Journal  d'Anal.  14  (1965),  161-169. 

9.  V.  I.  Arnold,  Sur  un  principe  variationnel  pour  lee  4cou laments  stationnaires  dee 
liquides  parfaits  et  ses  applications  aux  problames  de  stability  non  linrfares,  J. 
Mdcanique  5  (1966),  29-43. 

10.  V.  I.  Arnold,  Sur  la  gtometrie  differentia  He  dee  group  de  Lie  de  dismnsion  infinie  et 
ses  applications  a  l’hydrodynamique  das  fluids  parfait,  Ann.  Inst.  Fourier  (Grenoble) 
16  (1966),  319-361. 

11.  G.  B.  Whitham,  Variational  methods  and  applications  to  water  waves,  Proc.  Roy.  Soc. 
London,  A  299  (1967),  6-25. 

12.  J.  C.  Luke,  A  variational  principle  for  a  fluid  with  a  free  surface,  J.  Fluid  Mech.  17 
(1967),  395-397. 

13.  T.  B.  Benjamin,  Lectures  on  nonlinear  wave  motion,  Lectures  in  Applied  Mathematics, 
Vol.  15  (Ed.  A.  C.  Newell),  (User.  Math.  Society,  Providence,  1974. 


-46- 


»— arte  5.3.  The  resolution  of  problem  S  (1.32),  the  solitary  wave  problem,  is  carried 
out  by  letting  the  period  2k  approach  infinity  just  as  in  theorems  8.3  and  8.3 1  of  [1] 
for  4  •  0  or  4  >  0,  respectively.  The  assertions  of  theorem  1.2  concerning  exponential 
decay  and  the  case  k  ■  are  thus  covered,  completing  the  proof  of  the  theorem. 


-45- 


REFERENCES 


1.  R.  I.  l.  Turner,  Internal  waves  in  fluids  with  rapidly  varying  density,  Annali  della 
Scuola  Normals  -  Pisa,  Ser.  IV,  8  (1981),  513-573. 

2.  J.  L.  Bona,  D.  K.  Bose,  R.  E.  L.  Turner,  Finite  amplitude  steadywaves  In  stratified 
fluids,  Mathematics  Research  Center  TSR  Ho.  2401,  University  of  Wisconsin,  Madison, 
1981 i  Jour,  de  Math  Pure  at  Appl.,  to  appear. 

3.  J.  V.  Wahausen  and  S.  V.  Lai tone,  Surface  Waves,  Handbuch  der  Phyaik,  Vol.  IX, 

Strdmungsmechanik  III  (Ed.  S.  Flflgge)  Springer-Verlag,  Berlin,  1960. 

4.  C.  J.  Amick  and  J.  F.  Toland,  On  finite  amplitude  solitary  water  waves.  Arch.  Rat. 
Mech.  Anal.  76  (1980),  9-95. 

5.  A.  M.  Ter-Krikorov,  Thiorie  axacte  das  ondes  longues  statlonnaires  dans  un  liquids 

hdtSrogene  J.  Mtcanique  2  (1963),  351-376. 

6  K.  O.  Friedrichs  and  o.  H.  Hyers,  The  existence  of  solitary  waves.  Comm.  Pure  Appl. 
Math.  7  (1954),  517-550. 

7.  M.  A.  Lavrentiev,  A  contribution  to  the  theory  of  long  waves,  Amar.  Math.  Soc. 
Translation,  Ho.  102,  Providence,  RI,  1954. 

8.  P.  R.  Garabedian,  Surface  waves  of  finite  depth.  Journal  d’Anal.  14  (1965),  161-169. 

9.  V.  I.  Arnold,  Sur  un  principe  variatlonnel  pour  les  <cou laments  statlonnaires  des 
llquldes  parfaits  et  sas  applications  aux  problames  de  stability  non  lindares,  J. 
Mtcanique  5  (1966),  29-43. 

10.  V.  I.  Arnold,  Sur  la  gioaetrla  differential^  des  group  de  Lie  da  dimension  infinle  et 
ses  applicationa  a  l'hydrodynamique  des  fluids  parfait,  Ann.  Inst.  Fourier  (Grenoble) 
16  (1966),  319-361. 

11.  G.  B.  Whitham,  Variational  methods  and  applications  to  water  waves,  Proc.  Roy.  Soc. 
London,  A  299  (1967),  6-25. 

12.  J.  C.  Luke,  A  variational  principle  for  a  fluid  with  a  free  surface,  J,  Fluid  Mech.  17 
(1967),  395-397. 

13.  T.  B.  Benjamin,  Lectures  on  nonlinear  wave  motion.  Lectures  in  Applied  Mathematics, 
Vol.  15  (Ed.  A.  C.  Newell),  Amor.  Math.  Society,  Providence,  1974. 


-46- 


14.  T.  B.  Benjamin,  The  stability  of  solitary  waves,  Proc.  Roy.  Soc.  London,  A  328  (1972) 


153-183. 

15.  D.  D.  Bennett,  J.  L.  Bona,  R.  H.  Brown,  S.  E.  stansfield,  and  J.  D.  stronghair,  The 
stability  of  internal  solitary  waves  in  stratified  fluids  (to  appear). 

16.  J.  L.  Bona,  On  the  stability  theory  of  solitary  waves,  Proc.  Roy.  Soc.  London  A  344 
(1975),  363-374. 

17.  A.  Friedman  and  B.  Turkington,  Vortex  rings:  existence  and  asymptotic  estimates, 
Trans.  Amer.  Math.  Soc.  268  (1981),  1-37. 

18.  G.  Ready  and  J.  Norbury,  On  the  existence  theory  for  irrotational  water  waves.  Math. 
Proc.  Camb.  Phil.  Soc.  83  (1978),  137-157. 

19.  C.  J.  Amick,  L.  E.  Praenkel,  and  J.  P.  To land,  On  the  Stokes  conjecture  for  the  wave 
of  extreme  form,  Acta  Math.  148  (1982),  193-214. 

20.  M.  L.  Dubreil-Jacotin,  Sur  les  theoremes  d'existence  relatifs  aux  ondes  permanentes 
periodiques  a  deux  dimensions  dans  les  liquides  heterogenes,  J.  Math.  Pures  Appl.  19 
(1937),  43-67. 

21.  R.  R.  Long,  Some  aspects  of  the  flow  of  stratified  fluids.  Part  I.  A  theoretical 
investigation,  Tellus  5  (19S3),  42-57. 

22.  C.-S.  Yih,  On  the  flow  of  a  stratified  fluid,  Proc.  U.  S.  Nat.  Congr.  Appl.  Mech.  3rd 
(1958),  857-861. 

23.  T.  B.  Benjamin,  A  unified  theory  of  conjugate  flows,  Phil.  Trans.  Roy.  Soc.  London,  A 
269  (1971),  587-643. 

24.  0.  A.  Ladyzhenskaya  and  N.  N.  (Jral'tseva,  Linear  and  Quasilinear  elliptic  equations. 
Academic  Press,  New  York,  1968. 

25.  D.  Gilbarg  and  N.  S.  Trudinger,  Elliptic  Partial  Differential  Equations  of  Second 
Order,  Grundlehrer  der  Mathematische  wissenschaften  224,  Springer-Verlag,  Berlin, 
Heidelberg,  New  York,  1977. 

26.  J.  L.  Lions  and  I.  Magenes,  Problimes  aux  Limites  Non  Homogenes  et  Applications,  Vol. 


-47- 


1,  Dunod,  Paris,  1968 


27.  G.  Polya  and  G.  Sr *90,  Isoperimetric  Inequalities  in  Mathematical  Physics,  annals  of 


Mathematical  Studies  27,  Princeton  University  Press,  Princeton,  1950. 


RELT/ed 


REPORT  DOCUMENTATION  PAGE 


4.  TITLE  (mnd  Subtttlm) 

A  VARIATIONAL  APPROACH  TO  SURFACE  SOLITARY 
HAVES 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


1.  RECIPIENT'S  CATALOG  NUMBER 


S.  TYPE  OP  REPORT  *  PERIOD  COVEREO 

Summary  Report  -  no  specific 
_  reporting  period 


S.  PERFORMING  ORG.  REPORT  NUMBER 


7.  AUTHOR^ 

R.  E.  L.  Turner 


§.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Mathematics  Research  Center,  University  of 
610  Walnut  Street  Wisconsin 

Madison.  Wisconsin  53706 


<1.  CONTROLLING  OFFICE  NAME  AND  AOORESS 


i 


(See  Item  18  below) 


MONITORING  AGENCY  NAME  A  ADDRESSflf  tUllmrmnt  treat  Controlling  Olllcm) 


•  CONTRACT  OR  GRANT  NUMBERS; 

DAAG29-80-C-00  41 

MCS- 7904426 


10.  PROGRAM  element,  project,  task 

AREA  0  WORK  UNIT  NUMBERS 

Work  Unit  Number  1  - 
Applied  Analysis 


<2.  REPORT  DATE 

January  1983 


IS.  number  of  pages 
48 


IS.  SECURITY  CLASS,  (at  thin  import) 

UNCLASSIFIED 


«.  DISTRIBUTION  STATEMENT  (o I  Mia  Rmport) 

Approved  for  public  release;  distribution  unlimited. 


•7.  DISTRIBUTION  STATEMENT  (o I  aim  mbmtrmct  mntmrmd  In  Bloat  30.  II  Otltm rani  from  Rmpart) 


It.  SUPPLEMENTARY  NOTES 

U.  S.  Army  Research  Office 
P.  O.  Box  12211 
Research  Triangle  Park 
North  Carolina  27709 


It.  KEY  WORDS  (Co ntlnuo  on  eoroeoo  oldo  il  noeo 


National  Science  Foundation 
Washington,  DC  20550 


a  ary  mnd  Idmntlfy  by  black  numbmr) 


Solitary  wave,  surface  wave,  heterogeneous  fluid,  cnoidal  wave, 
critical  point,  symmetrization,  bifurcation 


20.  ABSTRACT  (Continue  an  raaaraa  alda  II  naeaaaary  and  Identity  by  block  numbmr) 

Two-dimensional  flow  of  an  incompressible,  inviscid  fluid  in  a  region 
with  a  horizontal  bottom  of  infinite  extent  and  a  free  upper  surface  is 
considered.  The  fluid  is  acted  on  by  gravity  and  has  a  non-diffusive, 
heterogeneous  density  which  may  be  discontinuous.  It  is  shown  that  the 
governing  equations  allow  both  periodic  and  single-crested  progressing  waves 
of  permanent  form,  the  analogues,  respectively,  of  the  classical  cnoidal  and 
solitary  waves.  The*-^  waves  are  shown  to  be  critical  points  of  flow  related 
functionals  and  are  reived  to  exist  by  means  of  a  variational  principle . 


EDITION  OF  I  NOV  ••  IB  OBBOLETE 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data  Bntered) 


to*  - 


